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1. Introduction 

Although the concept of Bose-Einstein |l|, ^ or intensity interferometry was dis- 
covered in particle and nuclear physics more than 30 years ago ||,^, some basic 
questions in the field are still unanswered, namely, what the form of the Bose- 
Einstein correlation functions is, and what this form means. However, even if the 
ultimate understanding of the effect is still lacking, the level of sophistication in the 
theoretical descriptions and the level of sophistication in the experimental studies 
of Bose-Einstein correlations and particle interferometry has increased drastically, 
particularly in the field of heavy ion physics |^ . 

1.1. W-mass determination and particle interferometry 

The study of Bose-Einstein correlations is interesting in its own right, but it should 
be noted that consequences may spill over into other fields of research, that are 
seemingly unrelated. Such is the topic of the W-mass determination at LEP2, a top 
priority research in high energy physics. It turned out that the non-perturbative 
Bose-Einstein correlations between the pions from decaying W+W~ pairs could be 
responsible for the presently largest systematic errors in W-mass determination at 
LEP2 |^,|^. Hence, the theoretical understanding and the experimental control of 
Bose-Einstein correlations at LEP2 is essential to make a precision measurement of 
the W mass, which in turn may carry information via radiative corrections about 
the value of the Higgs mass or signals of new physics beyond the Standard Model. 

1.2. Quark-gluon plasma and particle interferometry 

Heavy ion physics is the physics of colliding atomic nuclei. At the presently largest 
energies, the aim of heavy ion physics is to study the sub-nuclear degrees of free- 
dom by successfully creating and identifying the quark-gluon plasma (QGP). This 
presently only hypothetical phase of matter would consist of freely moving quarks 
and gluons, over a volume which is macroscopical relative to the characteristic 1 fm 
size of hadrons. 

Theoretically proposed signals of the expected phase transition from hot had- 
ronic matter to QGP were tested till now by fixed target experiments. At AGS, 
Brookhaven, collisions were made with nuclei as big as ^^^Au accelerated to 
14.5 AGeV bombarding energy. At CERN SPS, collisions were made with 60 and 
200 AGeV beams of nuclei, 200 AGeV beams of ^^S nuclei, 40 and 158 AGeV 
beams of ^osp]-, i^^clei ||^. The really heavy projectile runs were made relatively re- 
cently, the data are being published and the implications of the new measurements 
are explored theoretically, with claims of a possible QGP production at CERN SPS 
Pb -|- Pb reactions, however, without a clear-cut experimental proof of the identifi- 
cation of the new phase I). Both at CERN and at BNL, new collider experiments 
are planned and being constructed. The Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven will collide 100 -|- 100 AGeV ^^'^Au nuclei, which yields about 40 TeV 
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total energy in the center of mass frame. RHIC started to deliver its first results 
in 2000. The construction stage of the RHIC accelerator rings was declared to be 
complete by the US Department of Energy on August 14, 1999, during a NATO 
Advanced Study Institute in Nijmegen, The Netherlands, where the material of this 
review paper has been presented. The forthcoming Large Hadron Collider (LHC) 
at CERN is scheduled to start in 2005. LHC wiU coUide nuclei up to ^ospb with 
2.76 + 2.76 ATcV bombarding energy, yielding a total energy of 1150 TeV in the 
center of mass frame. The status quo has been summarized recently in Refs [^|-^|. 

At such large bombarding energies, the sub-nuclear structure of matter is ex- 
pected to determine the outcome of the experiments. However, the observed single 
particle spectra and two-particle correlations indicated rather simple dependences 
on the transverse mass of the produced particles [^;|l^, that had a natural ex- 
planation in terms of hydrodynamical parameterizations. Although hydrodynami- 
cal type of models are also able to fit the final hadronic abundances, spectra and 
correlations, [|| these models are not able to describe the ignition part of the pro- 
cess, thus their predictions are dependent on the assumed initial state. The hydro 
models come in two classes: i) hydro parameterizations, that attempt to parame- 
terize the flow, temperature and density distributions on or around the freeze-out 
hypersurfac e [|l6[-p4| by fitting the observed particle spectra and correlations, for 
example pj|, ^5| , |2§| -|3Cf|, but without solving the time-dependent (relativistic) hydro- 
dynamical equations. The class ii) comes in the form of hydrodynamical solutions, 
that assume an equation of state and an initial condition, and follow the time evolu- 
tion of the hydrodynamical system until a freeze-out hypersurface. These are better 
substantiated but more difficult to fit calculations, than class i) type of parameter- 
izations. The exact hydro solutions are obtained either in analytical forms, [^-Q, 
or from numerical solutions, see for example Refs pl]-^^. An even more substan- 
tiated approach is hydrodynamical approach with continuous emission of particles, 
which takes into account the small sizes of heavy ion reactions as compared to the 
mean free path of the particles ]45| ] . Such a continuous emission of hadrons during 
the time evolution of the hot and dense hadronic matter is supported by microscopic 
simulations |Q. 

In principle, the exact hydrodynamical solutions can be utilized in a time- 
reversed form: after fixing the parameters to describe the measured particle spectra 
and correlations at the time when the particles are produced, the hydro code can 
be followed backwards in time, and one may learn about the initial condition ||4^ 
in a given reaction: was it a QGP or a conventional hadron gas initial state? 



1.3. Basics of quantum statistical correlations 

Essentially, intensity correlations appear due to the Bose-Einstein or Fermi-Dirac 
symmetrization of the two-particle final states of identical bosons or fermions, in 
short, due to quantum statistics. 

The simplest derivation is as follows: suppose that a particle pair is observed, 
one with momentum fci the other with momentum fc2 ■ The amplitude has to be sym- 
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metrized over the unobservable variables, in particular over the points of emissions 
xi and X2. If Coulomb, strong or other final-state interactions can be neglected, 
the amplitude of such a final state is proportional to 

12 ^ L J, y J 

where + sign stands for bosons, — for fermions. If the particles are emitted in an 
incoherent manner, the observable two-particle spectrum is proportional to 

iV2(fci,fc2)cx J dxip{xi) j dX2p{x2) |Ai2|^ (2) 

and the resulting two-particle intensity correlation function is 

C2{kuk.) = , =i± \Pih - fc2)p (3) 

Ni[ki)N2[k2) 

that carries information about the Fourier-transformed space-time distribution of 
the particle emission 

p{q) ^ fdx e'«- pix) (4) 



as a function of the relative momentum q = ki — k2- 

As compared to the idealized case when quantum-statistical correlations are 
negligible (or neglected), Bose-Einstein or Fermi-Dirac correlations modify the 
momentum distribution of the hadron pairs in the final state by a weight factor 

(l±C0s[(ki-k2)-(xi-X2)]). 



1.4. Correlations between particle and heavy ion physics 

In case of pions, that are produced abundantly in relativistic heavy ion experi- 
ments, Bose-Einstein symmetrization results in an enhancement of correlations of 
pion pairs with small relative momentum, and the correlation function carries in- 
formation about the space-time distribution of pion production points. This in 
turn is expected to be sensitive to the formation of a transient quark-gluon plasma 
stage Q. 

In particle physics, reshuffling or modification of the momentum of pions in 
the fully hadronic decays of the W+W~ pairs happens due to the Bose-Einstein 
symmetrization of the full final stage, that includes symmetrization of pions with 
similar momentum from different W-s. As a consequence of this quantum inter- 
ference of pions, a systematic error as big as 100 MeV may be introduced to the 
W-mass determination from reconstruction of the invariant masses of (qq) systems 
in 4-jet events [^,0]. It is very difficult to handle the quantum interference of pions 
from the W"*" and jets with Monte-Carlo simulations, perturbative calculations 
and other conventional methods of high energy physics. 
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Unexpectedly, a number of recent experimental results arose suggesting that the 
Bose-Einstein correlations and the soft components of the single-particle spectra in 
high energy collisions of elementary particles show similar features to the same 
observables in high energy heavy ion physics |49|-^. 

These striking similarities of multi-dimensional Bose-Einstein correlations and 
particle spectra in high energy particle and heavy ion physics have no fully explored 
dynamical explanation yet. This review intends to give a brief introduction to 
various sub- fields of particle interferomctry, highlighting those phenomena that may 
have applications or analogies in various different type of reactions. The search for 
such analogies inspired a study of non-relativistic heavy ion reactions in the 30 - 
80 AMeV energy domain and a search for new exact analytic solutions of fireball 
hydrodynamics, reviewed briefly for a comparison. 

As some of the sections are mathematically more advanced, and other sections 
deal directly with data analysis, I attempted to formulate the various sections so 
that they be self-standing as much as possible, and be of interest for both the 
experimentally and the theoretically motivated readers. 



2. Formalism 



The basic properties of the Bose-Einstein n-particle correlation functions (BECF-s) 
can be summarized as follows, using only the generic aspects of their derivation. 
The n-particle Bose-Einstein correlation function is defined as 

where 7V„(ki, • • • ,k„) is the n-particle inclusive invariant momentum distribution, 
while 

1 d^"cr 

N^{^^,---,V^) = -E^,---E^- — (6) 

a dki • • • ak„ 

is the invariant n-particle inclusive momentum distribution. It is quite remarkable 
that the complicated object of Eq. (|^) carries quantum mechanical information on 
the phase-space distribution of particle production as well as on possible partial 
coherence of the source, can be expressed in a relatively simple, straight-forward 
manner both in the analytically solvable pion-laser model of Refs |52-56| as well as 
in the generic boosted-current formalism of Gyulassy, Padula and collaborators ||57| - 



as 



C„(ki,...,k„) = ^^::ip^ , (7) 

nG(k„ko 



i=l 



where tr*^"^ stands for the set of permutations of indices (1, 2, • • • , n) and Oi denotes 
the element replacing element i in a given permutation from the set of cr*^"' , and. 
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regardless of the details of the two different derivations, 



G(k„ k,) - ^Sk.i^k, (a^(k,)a(k,)) (8) 

stands for the expectation value of at(ki)a(kj). The operator a'l'(k) creates while 
operator a(k) annihilates a boson with momentum k. The quantity G'(ki, kj) corre- 
sponds to the first order correlation function in the terminology of quantum optics. 
In the boosted-current formalism, the derivation of Eq. (|^) is based on the assump- 
tions that i) the bosons are emitted from a semi-classical source, where currents are 
strong enough so that the recoils due to radiation can be neglected, ii) the source 
corresponds to an incoherent random ensemble of such currents, as given in a boost- 
invariant formulation in Ref. | |5^ , and Hi) that the particles propagate as free plane 
waves after their production. Possible correlated production of pairs of particles is 
neglected here. Note also the recent clarification of the proper normalization of the 
two-particle Bose-Einstein correlations [ ]60[ . 

A formally similar result is obtained when particle production happens in a 
correlated manner, generalizing the results of Refs p^-p6|,|6l|,|62|. Namely, the n- 
particle exclusive invariant momentum distributions of the pion-laser model read 
as 

n 

7V(")(ki, • • • , k„) = ^ n Gi(k„ k,J, (9) 

with 



Gi(k,k,) = ./Sk,£^k,Tr{piat(k,)a(kj)}, (10) 



where pi is the single-particle density matrix in the limit when higher-order Bose- 
Einstein correlations are negligible. Q.H. Zhang has shown |^^, that the ri-particle 
inclusive spectrum has a similar structure: 

n 

iV„(ki,...,k„) = ^ J|^(i^*'k„.)' (11) 

oo 

G(k„k,) = ^G„(k„k,). (12) 

n=l 

This result, valid only if the density of pions is below a critical value |^6|, was 
obtained if the multiplicity distribution was assumed to be a Poissonian one in the 
rare gas limit. The formula of Eq. (|lj) has been generalized by Q.H. Zhang in Ref. 
p3| to the case when the multiplicity distribution in the rare gas limit is arbitrary. 

The functions G„(ki, k^) can be considered as representatives of order n sym- 
metrization effects in exclusive events where the multiplicity is fixed to n, see Refs 
p3|-^ for more detailed definitions. The function G(ki,kj) can be considered 
as the expectation value of a^(ki)a(kj) in an inclusive sample of events, and this 
building block includes all the higher-order symmetrization effects. In the relativis- 
tic Wigncr function formalism, in the plane wave approximation G(ki,k2) can be 
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rewritten as 

G(ki,k2) = 5(1712,^^12) = J d'^x S{x,Ki2) exp{iqi2 ■ x) , (13) 

7^12 = 0.5(fci + A:2), (14) 
912 = h-h, (15) 

where a four- vector notation is introduced, k — (_Ek,k), and the energy of quanta 
with mass m is given by ~ y/ m? + , the mass-shell constraint. Notation 
a ■ b stands for the inner product of four-vectors. In the following, the relative 
momentum four- vector shall be denoted also as Afc = q = (go, fe, 9^, ^z) = (90, q), 
the invariant relative momentum is Q = ' 1- 

The covariant Wigner transform of the source density matrix, S{x,\ij is a 
quantum-mechanical analogue of the classical probability that a boson is produced 
at a given (x, fc) point in the phase-space, where x = {t,r) = (t,rx,ry,rz). The 
quantity S{x,Ki2) corresponds to the off-shell extrapolation of 5'(a::,k), as K^2 
y/ rn? + K^2- Fortunately, Bose-Einstein correlations are non-vanishing at small 
values of the relative momentum q, where — -E'ki2 ■ Due to the mass-shell 
constraints, G depends only on 6 independent momentum components. 

For the two-particle Bose-Einstein correlation function, Eqs ([7|,p|Jl3|) yield the 
following representation: 

C2(ki,k2) = l + i^^^i^i5i^. (16) 

Due to the unknown off-shell behavior of the Wigner functions, it is rather difficult 
to evaluate this quantity from first principles, in a general case. 

When comparing model results to data, two kind of simplifying approximations 
are frequently made: 

i) The on-shell approximation can be used for developing Bose-Einstein af- 
terburners to Monte-Carlo event generators, where only the on-shell part of the 
phase-space is modeled. In this approximation, Eq. ( p^ ) is evaluated with the on- 
shell mean momentum, K = {^J im? + K^j: K12). This on-shell approximation was 
used e.g. in Ref. |64| to sample S{x,K) from the single-particle phase-space distri- 
bution given by Monte-Carlo event generators, and to calculate the corresponding 
Bose-Einstein correlation functions in a numerically efficient manner. The method 
yields a straightforward technique for the inclusion of Coulomb and strong final- 
state interactions as well, see e.g. Ref. [Q. 

ii) The smoothness approximation can be used when describing Bose-Einstein 
correlations from a theoretically parameterized model, e.g. from a hydrodynamical 
calculation. In this case, the analytic continuation of S{x,'k) to the off-shell values 
of K is providing a value for the off-shell Wigner function S{x,Ki2). However, in 
the normalization of Eq. (^|), the product of two on-shell Wigner functions appear. 
In the smoothness approximation, one evaluates this product as a leading order 
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Taylor series in q of the exact expression S{0, K — q/2)S'(0, K + q/2). The resulting 
formula, 

relates the two-particle Bose-Einstein correlation function to the Fourier-trans- 
formed off-shell Wigner function S{x,K). This provides an efficient analytic or 
numeric method to calculate the BECF from sources with known functional forms. 
The correction terms to the smoothness approximation of Eq. are given in Ref. 
[ p3| . These corrections are generally on the 5% level for thermal like momentum 
distributions. 



3. Model-Independent Analysis of Short-Range Correlations 

Can one model-independently characterize the shape of two-particle correlation func- 
tions? Let us attempt to answer this question on the level of statistical analysis, 
without theoretical assumptions on the thermal or non-thermal nature of the par- 
ticle emitting source. In this approach, the usual theoretical assumptions are not 
made, neither on the presence or the negligibility of Coulomb and other final-state 
interactions, nor on the presence or the negligibility of a coherent component in 
the source, nor on the presence or the negligibility of higher-order quantum statis- 
tical symmetrization effects, nor on the presence or the negligibility of dynamical 
effects (e.g. fractal structure of gluon-jets) on the short-range part of the correlation 
functions. The presentation follows the lines of Ref. ||6^. The reviewed method is 
really model-independent, and it can be applied not only to Bose-Einstein correla- 
tion functions but to every experimentally determined function, which features the 
properties i) and ii) listed below. 

The following experimental properties are assumed: 

i) The measured function tends to a constant for large values of the relative 
momentum. 

ii) The measured function has a non-trivial structure at a certain value of its 
argument. 

The location of the non-trivial structure in the correlation function is assumed 
for simplicity to be close to Q = 0. 

The properties i) and ii) are well satisfied by e.g. the conventionally used two- 
particle Bose-Einstein correlation functions. For a critical review on the non-ideal 
features of short-range correlations, (e.g. non-Gaussian shapes in multi-dimensional 
Bose-Einstein correlation studies), we recommend Ref. [ |66[ . 

The core/halo intercept parameter A* is defined as the extrapolated value of 
the two-particle correlation function at Q — 0, see Section || for greater details. 
It turns out that A^, is an important physical observable, related to the degree of 
partial restoration of Ua{^) symmetry in hot and dense hadronic matter [p^,p8|, as 



reviewed in Section 15 
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Various non-ideal effects due to detector resolution, binning, particle mis-identi- 
fication, resonance decays, details of the Coulomb and strong final-state interactions 
etc. may influence this parameter of the fit. One should also mention, that if all 
of these difficulties are corrected for by the experiment, the extrapolated intercept 
parameter A, for like-sign charged bosons is (usually) not larger, than unity as a 
consequence of quantum statistics for chaotic sources, even with a possible admix- 
ture of a coherent component. However, final-state interactions, fractal branch- 
ing processes of gluon jets, or the appearance of one-mode or two-mode squeezed 
states |^,|7^ in the particle emitting source might provide arbitrarily large values 
for the intercept parameter. 

A really model-independent approach is to expand the measured correlation 
functions in an abstract Hilbert space of functions, ft is reasonable to formulate 
such an expansion so that already the first term in the series be as close to the 
measured data points as possible. This can be achieved if one identifies |^,|7l[] the 
approximate shape (e.g. the approximate Gaussian or the exponential shape) of the 
correlation function with the abstract measure fi{t)dt in the abstract Hilbert-space 
Ti.. The orthonormality of the basis functions 0„ (t) in Ti can be utilized to guarantee 
the convergence of these kind of expansions, see Refs [^,0 for greater details. 

3.1. Laguerre expansion and exponential shapes 

If in a zeroth order approximation the correlation function has an exponential shape, 
then it is an efficient method to apply the Laguerre expansion, as a special case of 
the general formulation of Refs |6^, Q : 



In this and the next subsection, Q stands symbolically for any, experimentally 
chosen, one dimensional relative momentum variable. The fit parameters are the 
scale parameters Af, Xl, Rl and the expansion coefficients ci, C2, ... . The nth 
order Laguerre polynomials are defined as 



C2{Q) = AA{ 



I + Xl exp(-Qi?i) [l + ciLi{QRl) -f ^L2{QRl) + ..•] } 



(18) 



L„(i) = exp(i)^t"exp(-i), 
they form a complete orthogonal basis for an exponential measure as 



(19) 



<5, 



'n.m 



OC 




(20) 



The first few Laguerre polynomials are explicitly given as 



Lo{t) = 1, 

Li(t) = t-1, 

L2{t) = i2_4t + 2, .... 



(21) 
(22) 
(23) 
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As the Laguerre polynomials are non- vanishing at the origin, C{Q = 0) =/= l+A^. 
The physically significant core/halo intercept parameter A* can be obtained from 
the parameter Al of the Laguerre expansion as 

A, = Ai[l-ci+c2 -...]. (24) 



3.2. Edgeworth expansion and Gaussian shapes 



If, in a zeroth-order approximation, the correlation function has a Gaussian shape, 
then the general form given in Ref. |^2| takes the particular form of the Edgeworth 
expansion fTlnZl] 



C(Q) 



7V{1 + Ab exp{~Q^Rl 



1 + -±H^{V2QRe) + -^Hi{V2QRE) 



(25) 



The fit parameters are the scale parameters A/", Xe, Re, and the expansion coeffi- 
cients K3, K4, ... that coincide with the cumulants of rank 3, 4, ... of the correlation 
function. The Hermite polynomials are defined as 



Hn{t) 



exp(t2/2) 



d 
dt 



exp(-tV2), 



(26) 



they form a complete orthogonal basis for a Gaussian measure as 

/oo 
dt exp{-t^/2)H,,{t)H„,{t). (27) 
-OO 



The first few Hermite polynomials are listed as 



Hiit) 
H2(t) 
Hz(t) 



t" 



1, 

it, 

6t' 



(28) 
(29) 
(30) 
(31) 



The physically significant core/halo intercept parameter A* can be obtained from 
the Edgeworth fit of Eq. (H) as 



A* 



A, 



(32) 



This expansion technique was applied in the conference contributions ||7l|,|72| to the 
AFS minimum bias and 2-jet events to characterize successfully the deviation of 
data from a Gaussian shape. It was also successfully applied to characterize the 
non-Gaussian nature of the correlation function in two-dimensions in case of the 
preliminary E802 data in Ref. |7lf|, and it was recently applied to characterize the 
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non-Gaussian nature of the three-dimensional two-pion BECF in e++ e reactions 
at LEPl 

Figure 111 indicates the ability of the Laguerre expansions to characterize two 
well-known, non-Gaussian correlation functions Jsst : the second-order short-range 
correlation function _D| (Q) as determined by the UAl and the NA22 experiments |7^, 
[TsI . The convergence criteria of the Laguerre and the Edgeworth expansions is given 
in Ref. |6|. 



Laguerre expansion fit 




10"^ 10^ 10^ 1 10"^ 10^ 10^ 1 

UAl data NA22 data 



Fig. 1. Laguerre expansion of NA22 and UAl short-range correlations is shown 
by the solid line. Dashed line stands for the best exponential fit, which clearly 
underestimates the strength of the measured points at low values of the squared 
invariant momentum difference = — (fci — fc2)^- (Note the logarithmic horizontal 
and vertical scales.) 

From Table ^ the core/halo model intercept parameter is obtained as A* = 
1.14 ±0.10 (UAl) and A* = 1.11 ±0.17 (NA22). As both of these values are within 
errors equal to unity, the maximum of the possible value of the intercept parameter 
A* in a fully chaotic source, we conclude that either there are other than Bose- 
Einstein short-range correlations observed by both collaborations, or the full halo 
of long lived resonances is resolved in case of this measurement [7^-[79| . 

If the two-particle BECF can be factorized as a product of (two or more) func- 
tions of one variable each, then the Laguerre and the Edgeworth expansions can be 
applied to the multiplicative factors — functions of one variable, each. This method 
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Table 1. Laguerre fits to UAl and NA22 two-particle correlations 





UAl 


NA22 




1.355 ± 0.003 


0.95 ± 0.01 


\l 


1.23 ± 0.07 


1.37 ± 0.10 


Rl [fm] 


2.44 ± 0.12 


1.35 ± 0.14 


Cl 


0.52 ± 0.03 


0.63 ± 0.06 




0.45 ± 0.04 


0.44 ± 0.06 


X^/NDF 


41.2/41 = 1.01 


20.0/34 = 0.59 



was applied recently to study the non-Gaussian features of multi-dimensional Bose- 



Einstein correlation functions e.g. in Refs |52, 72[. The full, non-factorized form of 
two-dimensional Edgeworth expansion and the interpretation of its parameters is 
described in the handbook on mathematical statistics by Kendall and Stuart |Q . 

4. Coulomb Wave Corrections for Higher-Order Correlations 

The short-range part of the two- and multi-particle correlation function of charged 
particles is strongly effected by Coulomb interactions. Even in the non-rclativistic 
case, the n-body Coulomb scattering problem is solvable exactly only for the n = 2 
case, the full 3-body Coulomb wave-function is unknown. However, when studying 
higher-order Bose-Einstein correlations and e.g. searching for the onset of (partial) 
coherence in the source, it is desired that the Coulomb-induced correlations be 
removed from the data. 

In any given frame, the boost-invariant decomposition of Eq. ([l3| ) can be rewrit- 
ten into the following, seemingly not invariant form: 

G(ki,k2) = y"d3x ^Ki^W exp(iqi2x), (33) 

^Ki.(x) = y"dt exp(-^/3,,^^ql2^) 5(x,i,ifi2), (34) 
/3^^^ = (ki+k2)/(i?i+i?2). (35) 

Note that the relative source function 5ki2(x) reduces to a simple time integral 
over the source function S{x, K) in the frame where the mean momentum of the 
pair (hence the pair velocity Pki^) vanishes. 

Based on a Poisson cluster picture, the effect of multi-particle Coulomb final- 
state interactions on higher-order intensity correlations is determined in general in 
Ref. with the help of a scattering wave function which is a solution of the 
n-body Coulomb Schrodinger equation in (a large part of) the asymptotic region of 
the n-body configuration space. 

If n particles are emitted with similar momenta, so that their n-particle Bose- 
Einstein correlation functions may be non-trivial, Eqs (^3|-^5|) form the basis for 
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Fig. 2. Coulomb wave- function correction factor and generalized Gamow correction 
factor for 5-particle correlation functions, for a Gaussian source with Rq = 5 fm 

evaluation of the Coulomb and strong final-state interaction effects on the observ- 
ables for any given cluster of particles, assuming that the relative motion of the 
particles is non-relativistic within the cluster, see Ref. Q. The Coulomb correc- 
tion factor can be integrated for arbitrary large number of particles and for 
any kind of model source, by replacing the plane wave approximation with the ap- 
proximate n-body Coulomb wave- function. In the limit of vanishing source sizes, 
the generalization of the Gamow penetration factor was obtained to the correlation 
function of arbitrary large number of particles In particular, Coulomb effects 
on the n-particle Bose-Einstein correlation functions of similarly charged particles 
were studied for Gaussian effective sources, for n = 3 in Ref. |^ and for n = 4 
and 5 in Ref. . For the typical i? = 1 fm effective source sizes of the elementary 
particle reactions, the generalized n-body Gamow penetration factor gave rather 
precise estimates of the Coulomb correction (within 5% from the Coulomb-wave 
correction). In contrast, for typical effective source sizes observed in high energy 
heavy ion reactions. Fig. |^ indicates that the new Coulomb wave-function inte- 
gration method allows for a removal of a systematic error as big as 100% from 
higher-order multi-particle Bose-Einstein correlation functions. See Ref. |^l| for 
greater details. 
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5. Core/Halo Picture of Bose— Einstein Correlations 

The core/halo model [|76|,|83|-|8^ deals with the consequences of a phenomenological 
situation, when the boson source can be considered to be a superposition of a 
central core surrounded by an extended halo. In the forthcoming sections, final- 
state interactions are neglected, we assume that the data are corrected for final-state 
Coulomb (and possibly strong) interactions. 

Bose-Einstein correlations are measured at small relative momenta of particle 
pairs. In order to reliably separate the near-by tracks of particle pairs in the region of 
the Bose enhancement, each experiment imposes a cut-off Qmin, the minimum value 
of the resolvable relative momentum. The value of this cut-off may vary slightly 
from experiment to experiment, but such a cut-off exists in each measurement. 

In the core/halo model, the following assumptions are made: 

Assumption 0: The emission function does not have a no-scale, power-law-like 
structure. This possibility was discussed and related to intermittency and effective 
power-law shapes of the two-particle Bose-Einstein correlation functions in Ref. . 

Assumption 1: The bosons are emitted either from a central part or from the 
surrounding halo. Their emission functions are indicated by S'c(x,k) and S';i(x,k), 
respectively. According to this assumption, the complete emission function can be 
written as 



S{x,k)^Sc{x,k) + Sh{x,k), (36) 

and Slx.k) is normalized to the mean multiplicity, J d*x {dk/ E)S{x,'k) = (n). 

Assumption 2: The emission function that characterizes the halo is assumed to 
change on a scale Rh that is larger than i^max ~ ^/Qmim the maximum length-scale 
resolvable by the intensity interferometry microscope. The smaller central core 
of size Rc is assumed to be resolvable, 

Rh > Rinax > Rc ■ (37) 

This inequality is assumed to be satisfied by all characteristic scales in the halo and 
in the central part, e.g. in case the side, out or longitudinal components [^,^ of 
the correlation function are not identical. 

Assumption 3: The core fraction /c(k) — Ncik)/Ni{k) varies slowly on the 
relative momentum scale given by the correlator of the core |8^] . 

The emission function of the core and the halo are normalized as 

j d^x^S,{x,k) = (n), and j d^x^Sh{x,k) = (n)^. (38) 
One finds [fzl^ that 



A^i(k) = iV,(k)+A^,,(k) and {n) = {n), + {n)h. 



(39) 
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Q /MeV/ 

Fig. 3. The shape of the BECF is illustrated for a source containing a core 
and a large halo. The contribution from the halo is restricted to the shaded area, 
while the shape of the BECF outside this interval is determined completely by the 
contribution of the core. If the resolution for a given experiment is restricted to 
Q > 10 MeV, then an effective and momentum dependent intercept parameter, 
A*(y,rnt) will be measured, which can be combined with the measured momentum 
distribution to determine the the momentum distribution of the particles emitted 
directly from the core. 

Note that in principle the core as well as the halo part of the emission function 
could be decomposed into more detailed contributions, e.g. 

5,(a;,k) = J2 ^hl(^,k). (40) 

r— a;, 77,77', A"^ 

In case of pions and NA44 acceptance, the lu mesons were shown to contribute to 
the halo, Ref. |7^. For the present considerations, this separation is indifferent, 
as the halo is defined with respect to Qmim the experimental two-track resolution. 
For example, if Qmin = 10—15 MeV, the decay products of the cu resonances 
can be taken as parts of the halo |7|]. If future experimental resolution decreases 
below 5 MeV and the error bars on the measurable part of the correlation function 
decrease significantly in the Q < fi/Vuj = 8 MeV region, the decay products of the lo 
resonances will contribute to the resolvable core, see Refs [Q|7^ for greater details. 

If Assumption 3 is also satisfied by some experimental data set, then Eq. ( [T^ ) 
yields a particularly simple form of the two-particle Bose-Einstein correlation func- 
tion: 
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C.(k„k.) = 1 + A4K) (41) 

S^c(0,ki)5c(0,k2) 



i+A.(K) ';;^r:;.^' > (42) 



^c(Afc,g) 

|5e(0,if)|^ 

where mean and the relative momentum four-vectors are defined as 

A' = 0.5(fci + fc2), Afc = fci-fc2, (43) 

with K = (X^jK) and Ak — (A/c^jAk), and the effective intercept parameter 
A,(K) is given as 

A,(K) = [7V,(K)/iVl(K)]^ (44) 

As emphasized in Ref. [ f76{ , this effective intercept parameter A* shah in general 
depend on the mean momentum of the observed boson pair, which within the er- 
rors of Qmin coincides with any of the on-sheU four-momentum ki or k2 . Note that 
A* ^ Xxct — 1, the latter being the exact intercept parameter at Q = MeV. 
The core/halo model is summarized in Fig. ^, see Ref. ||7^ for further details. 
The core/halo model correlation function is compared to the so-called "model- 
independent" , Gaussian approximation of Refs and to the full correlation 
function in Fig. ^, see appendix of Ref. and that of Ref. Q for further details. 

The measured two-particle BECF is determined for | Ak |> Qmin ~ 10 MeV/c, 
and any structure within the | Ak |< Qmin region is not resolved. However, the 
(c, h) and {h, h) type boson pairs create a narrow peak in the BECF exactly in this 
Ak region according to Eq. (p6|), which cannot be resolved according to Assump- 
tion 2. 

The general form of the BECF of systems with large halo, Eq. (^2|), coincides 
with the most frequently applied phenomenological parameterizations of the BECF 
in high energy heavy ion as well as in high energy particle reactions ]89| ]. Previ- 
ously, this form has received a lot of criticism from the theoretical side, claiming 
that it is in disagreement with quantum statistics [ pO[ or that the A parameter is 
just a kind of fudge parameter, "a measure of our ignorance". In the core/halo 
picture, Eq. (^) is derived with a standard inclusion of quantum statistical effects. 
Reactions including e^-f e^ annihilations, lepton-hadron and hadron-hadron reac- 
tions, nucleon-nucleus and nucleus-nucleus collisions are phenomenologically well 
described [p^ by a core/halo picture. 



5.1. Partial coherence and higher-order correlations 

In earlier studies of the core/halo model [^|8^ it was assumed that Sc{x,p) de- 
scribes a fully incoherent (thermal) source. In Ref. an additional assumption 
was also made: 

Assumption 4- A part of the core may emit bosons in a coherent manner: 



Sc{x,k)=SP{x,k) + Sl{x,k), 



(45) 
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Fig. 4. Comparison of the full correlation function (full line) to the core/halo model 
approximation (dashed line) and to the "model-independent" Gaussian approxima- 
tion (dotted line) 

where upper index p stands for coherent component (which leads to partial coher- 
ence), upper index i stands for incoherent component of the source. 
The invariant spectrum is given by 



One can introduce the momentum dependent core fractions /c(k) and partially 
coherent fractions p(k) as 




(46) 



and the core contribution is a sum: 




(47) 



/e(k) = Ar,(k)/Ar(k), 
Pc(k) = NP{k)/Nc{k). 



(48) 
(49) 



Hence the halo and the incoherent fractions fh,Pi are 



A(k) = iV;,(k)/iV(k) = 1 - /e(k), 
/,(k) = 7V^(k)/Are(k) = l-pe(k). 



(50) 
(51) 
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5.2. Strength of the n-particle correlations 

We denote the n-particle correlation function of Eq. (||) as 

C„(l,2,...,n) . a.(k.k.,...,kj.^^^^;^^^^ (52) 

where a symbolic notation for is introduced, only the index of k is written out in 
the argument. In the forthcoming, we shall apply this notation consistently for the 
arguments of various functions of the momenta, i.e. /(k^, k^, k^) is symbolically 
denoted by /(i, j, to). 

The strength of the n-particle correlation function (extrapolated from a finite 
resolution measurement to zero relative momentum for each pair) is denoted by 
C„(0), given 6^ by the following simple formula, 



(53) 



Here, aj indicates the number of permutations, that completely mix exactly j non- 
identical elements. There are (") different ways to choose j different elements from 
among n different elements. Since all the n! permutations can be written as a sum 
over the fully mixing permutations, the counting rule yields a recurrence relation 
for aj, Refs p5|,p6t: 

an = ^--^i (54) 
ao = 1. (55) 



The first few values of Uj are given as 

ai = 0, a2 — 1, as = 2, a4 = 9, — 44, ag — 265, (56) 
the first few intercept parameters, A,.„ — C„(0) — 1, are given as 

A,,2 = /2[(l-pc)' + 2pc(l-Pc)], (57) 

+2f^[il-Pcf + Spcil^Pc)% (58) 

A,.4 = 6f^[{l-pcf + 2p,{l-p,)] 

+8f![{l^Pcf + 3pc{l~Pc)^] 

+9f^[il-p,f + Ap,{l~p,)% (59) 

A*.5 = 10/2[(l-p,)2 + 2p,(l-p,)] 

+20/3[(l-p,)3 + 3pc(l-Pc)'] 
+45/4[(l-p,)4 + 4p,(l-p,)3] 

+Uf^[il-P,f + 5pc{l-pc)% (60) 
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In general, terms proportional to in the incoherent case shall pick up an 
additional factor [(1 — PcY + jPc{^ — Pc)"'^^] in case the core has a coherent com- 
ponent |]85|,|8^. This extra factor means that either all j particles must come from 
the incoherent part of the core, or one of them must come from the coherent, the 
remaining J — 1 particles from the incoherent part. If two or more particles come 
from the coherent component of the core, the contribution to intensity correlations 
vanishes as the intensity correlator for two coherent particles is zero psf . 

If the coherent component is present, one can introduce the normalized incoherent 
and partially coherent core fractions as 

Slij^k) = (61) 

~4U,k) = f^. (62) 

In the partially coherent core/halo picture, one obtains the following closed form 
for the order n Bose-Einstein correlation functions : 

n n , ( j 

C„(l,...,n) = 1 + ^ ^ ^{\Yfc{mk)[l-pc{mk)]sl{mk,rnp^) 

j=2 mi ...mj = l p(j) 

j 3 \ 

+ '^fc{m)Pc{mi) sP{mi,mp^) Y\ fc{mk)[l ~ Pc{mk)] sl{mk,mpj > . (63) 

/=! k=l,k^l J 

Here, p^^^ stands for the set of permutations that completely mix exactly j ele- 
ments. Pi stands for the permuted value of index i in one of these permutations. 
By definition, pi ^ i for all i = 1,2,..., j. The notation S' indicates summation 
for different values of indexes, nii ^ mi for all i, I pairs. The expression Eq. ( |6^ ) 
contains two (momentum dependent) phases in the Fourier-transformed, normal- 
ized source distributions: one denoted by 0*(k,„,k„) in the Fourier-transformed 
normalized incoherent core emission function, s^(km,k„) and another independent 
phase denoted by 0'^(km,k„) is present in the the Fourier-transformed normalized 
coherent core emission function, s^(km,k„). One can write 

SJ,(k™,k„) = |S^(k„,k„)|exp[i0''(k,„,k„)], (64) 
5?(k„,k„) = |s?(k„,k„)|exp[i0''(k„,k„)]. (65) 

The shape of both the coherent and the incoherent components is arbitrary, but 
corresponds to the space-time distribution of particle production. If the variances 
of the core are finite, the emission functions can be parameterized by Gaussians, for 
the sake of simplicity |78). If the core distributions have power-law-like tails, like 
in case of the Lorentzian distribution [flSl, then the Fourier-transformed emission 



20 



T. Csorgo 



functions correspond to exponentials or to power-law structures. For completeness, 
we list these possibilities below: 

(66) 
(67) 
(68) 





= exp(-i?2Q^^) 


or 


ScO^rm k„)| 


= exp(-i?iQ,„„) 


or 


Sc(km, k,j)| 




etc 


ScO^rm kri)| 




or 


^c(^™i kn)| 


= exp(-i?p(5,„„) 


or 


Sc(kmi k„)p 




etc 



(69) 
(70) 
(71) 

In the above equations, subscripts i and p index the parameters belonging to the 
incoherent or to the partially coherent components of the core, and Qmn stands for 
certain experimentally defined relative momentum component determined from km 
and k„. 

A straightforward counting yields that in the limiting case when all momenta 
are equal, the simple formula of Eq. (|53) follows from the shape of the n-particle 



Bose-Einstein correlation functions of Eq. (63), as s\{i,i) — sP{i,i) — 1 



5.3. Graph rules 

Graph rules were derived for the evaluation of the n-particle correlation function 
C„(ki, k„) in Ref. Graphs contributing to the n ~ 2 and 3 case are shown 
in Fig. ^, the case of n = 4 is shown in Fig. 

Circles can be either open or full. Each circle carries one label (e.g. j) standing 
for a particle with momentum k^ . Full circles represent the incoherent core compo- 
nent by a factor /c(i)[l — Pc(j)], whereas open circles correspond to the coherent 
component of the core, a factor of fc(j)pc(j)- 

For the n-particle correlation function, all possible j-tuples of particles have to 
be found. Such j-tuples can be chosen in (") different manner. In a j-tuple, either 
each circle is filled, or the circle with index k is open and the other j — \ circle 
is filled, which gives j + 1 different possibilities. All the permutations that fully 
mix either j = 2 or 3, or n different elements have to be taken into account for 
each choice of filling the circles. The number of different fully mixing permutations 
that permute the elements ii, is given by aj and can be determined from the 
recurrence of Eq. (|^ . 

Lines, that connect a pair of circles (or vertexes) stand for factors that 

depend both on k^ and kj. Full lines represent incoherent-incoherent particle pairs, 
and corresponds to a factor of s],{i,j). Dashed lines correspond to incoherent- 
coherent pairs, and carry a factor of s^(j,j). The lines are oriented, they point 
from circle i to circle j, corresponding to the given permutation, that replaces 
element j by element i. Dashed lines start from an open circle and point to a full 
circle. 
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• = f,(k,)[l-p^(k,)] i j = sUk„kj) 

o = f^(k,)pi(ki) i--^- j =sj(ki,kj) 



I I 




Fig. 5. Graphs determining the second and the third order correlation function for 
partially coherent core/halo sources 

All graphs contribute to the order n correlation function, that are in agreement 
with the above rules. The result corresponds to the fully mixing permutations 
of all possible j'-tuples (j = 2, ...,n) chosen in all possible manner from elements 
(l,2,...,n). 

Each graph adds one term to the correlation function, given by the product 
of all the factors represented by the circles and lines of the graph. Note that the 
directions of the arrows matter. The correlation function C(l, ...,n) is given by 1 
plus the sum of all the graphs. 

Note that for the n-particle cumulant correlation function, n circles, represent- 
ing the n particles, should be connected in all possible manner corresponding only 
to the fully mixing permutations of elements (1, ...,n). Disconnected graphs do not 
contribute to the cumulant correlation functions, as they correspond to permuta- 
tions, that either do not mix all of the n elements or can be built up from two or 
more independent permutations of certain sub-samples of elements (1,2, n). 

5.4. Application to three-particle correlation data 

In the CERN SPS S -I- Pb reactions, the strength of the two- and three-particle 
correlation functions was determined experimentally by the NA44 Collaboration as 
K.2 = 0.44 ± 0.04 in Ref. |l| and by A,,3 = 1-35 ± 0.12 ± 0.09, Ref. 0. Note 
that the value of A*_3 was determined with the help of the Coulomb 3-particle wave- 
function integration method of Ref. Q , reviewed in Section |[ because the estimate 
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Fig. 6. Graphs determining the fourth order correlation function for partially 
coherent core/ halo sources 
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Fig. 7. Allowed regions for possible values of the core fraction fc and the partially 
coherent fraction pc are evaluated on the two standard deviation level from the 
intercept parameter of the second and the third order BE correlation functions, 

A2,* and As^, 

based only on the 3-body Gamow penetration factor introduced unacceptably large 
systematic errors to the three-particle Bose-Einstein correlation function. 

The two experimental values, A*, 2 and A*. 3 can be fitted with the two theoretical 
parameters fc and Pc, as done in Ref. Figure [t] illustrates the 2 a contour plots 
in the {fc,Pc) plane, obtained using the published value of A, .2 — 0.44 ± 0.04 and 
the preliminary value of A, ^3 = 1.35 ± 0.12. A range of (fcPc) values is found 
to describe simultaneously the strength of the two-particle and the three-particle 
correlation functions within two standard deviations from these values. Thus neither 
the fully chaotic, nor the partially coherent source picture can be excluded at this 
level of precision. 

Table 2. Strength of higher-order correlation functions for various core fractions 
and partially coherent fractions allowed by NA44 2- and 3-particle correlation data 



fc 


Pc 


A*^2 


A*, 3 


A*, 4 


A*. 5 


0.60 


0.00 


0.36 


1.51 


5.05 


17.17 


0.70 


0.50 


0.37 


1.45 


4.25 


11.87 


1.00 


0.75 


0.44 


1.63 


4.33 


10.47 



Cramer and Kadija pointed out, that for higher values of n the difference be- 
tween a partially coherent source and between the fully incoherent particle source 
with an unresolvable component (halo or mis-identified particles) will become larger 
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and larger Indeed, similar values can be obtained for the strength of the second 
and third order correlation function, if the source is assumed to be fully incoherent 
(/c = 0.6, Pc = 0) or if the source has no halo but a partially coherent component 
(/c = l,Pc = 0.75), but the strength of the 5th order correlation function is almost 
a factor of 2 larger in the former case, as can be seen from Table |^. Precision mea- 
surements of 4th and 5th order correlations are necessary to determine the value of 
the degree of partial coherence in the pion source. 

6. Particle Interferometry in e++ Reactions 

The hadronic production in e"'"e~ annihilations is usually considered to be a basi- 
cally coherent process and therefore no Bose-Einstein effect was expected, whereas 
hadronic reactions should be of a more chaotic nature giving rise to a sizable effect. 
It was even argued that the strong ordering in rapidity, preventing neighboring 
TT~TT~ or 7r+7r+ pairs, would drastically reduce the effect j9^. Therefore it was a 
surprise when G. Goldhaber at the Lisbon Conference in 1981 |94j presented data 
which showed that correlations between identical particles in e'^e^ annihilations 
were very similar in size and shape to those seen in hadronic reactions, see the 
review paper Ref. |^ for further details. 

6.1. The Andersson-Hofmann model 

The Bose-Einstein correlation effect, a priori unexpected for a coherent process, 
has been given an explanation within the Lund string model by B. Andersson and 
W. Hofmann The space-time structure of an e''"e~ annihilation is shown for 
the Lund string model |9^] in Fig. ||. The probability for a particular final state is 
given by the expression 

Prob. ~ phasespace • exp(— (72) 

where A is the space-time area spanned by the string before it breaks and & is a 
parameter. The classical string action is given by 5 = kA, where k is the string ten- 
sion. It is natural to interpret the result in Eq. ( [72|) as resulting from an imaginary 
part of the action such that 

S={k + lb/2) a, (73) 

and an amplitude M given by 

M - exp(iS'), (74) 

which implies 

Prob. -I M exp(-M). (75) 

Final states with two identical particles are indistinguishable and can be ob- 
tained in different ways. Suppose that the two particles indicated as 1 and 2 in 
Fig. ^ are identical, then the hadron state in the left panel can be considered as 
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Fig. 8. Andersson-Hofmann interpretation of Bose-Einstein correlations in the 
Lund string model. Ai^2 denotes the space-time area of a color field enclosed by 
the quark loop in c+e~ annihilation. Two particles 1 and 2 are separated by the 
intermediate system /. When the particles 1 and 2 are identical, the configuration 
in the left side is indistinguishable from that of the right side, and their amplitudes 
for production must be added. The probability of production will depend on the 
difference in area = Ai — A2, shown as the hatched area. 



being the same as that in the right panel (where 1 and 2 arc interchanged). The 
amplitude should, for bosons, be the sum of two terms 

M ~ exp[i(K + ib/2)Ai] + exp[i(K + ib/2)A2] , (76) 

where Ai and A2 are the two string areas, giving a probability proportional to 

I M [exp(-6Ai) + exp{~bA2)] ■ 

with A A = Ai — A2. The magnitudes of k and b arc known from phenomenological 
studies. The energy per unit length of the string is given by k ~ 1 GeV/fm, 
and b describes the breaking of the string at a constant rate per unit area, b/K^ « 
0.7 GeV~^ 1^^. The difference in space-time area A^l is marked as the hatched area 
in Fig. ||. It can be expressed by the (t, r^) components {E, k) of the four- momenta 
of the two identical particles 1 and 2, and the intermediate system /: 

= [E2ki - Eik2 + Ej{ki ~ fca) - kj{Ei - E2)]/k^ . (78) 

To take into account also the component transverse to the string a small addi- 
tional term is needed. The change in area A^ is Lorentz invariant to boosts 
along the string direction and is furthermore approximately proportional to Q = 

V-(fci-fe)'- 

The interference pattern between the amplitudes will be dominated by the 
phase change of A$ — kAA. It leads to a Bose-Einstein correlation which, as 
a function of the four-momentum transfer, reproduces the data well but shows a 



cos{kAA) 
cosh(6AA/2) 



(77) 
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steeper dependence at small Q than a Gaussian function. A comparison to TPC 
data confirmed the existence of such a steeper than Gaussian dependence on Q, 
although the statistics at the small Q- values did not allow a firm conclusion |]89|,p7|. 

Recently, the interest for multi-dimensional analysis of Bose-Einstein correla- 
tions increased also in the particle physics community, see Ref. for a critical 
review of the present status. 

I would like to highlight three interesting features: i) The effect seems to de- 
pend on the transverse momentum of the produced pion pairs, decreasing effective 
radii were observed for increasing transverse mass pi , 52 1 . This effect is also seen 



in the LUBOEI algorithm of JETSET, although no intrinsic momentum-dependent 
scale is plugged into the algorithm [Q. ii) The three-dimensional Bose-Einstein 
correlations of L3 indicate a non-Gaussian structure |52|. in) The effective source 
sizes of heavier particles (K, A) were measured recently [Q, based on spin statis- 
tics developed by Alexander and Lipkin |100|. The measured source sizes show a 
clear decrease with increasing particle masses. The latter effect was explained by 
Alexander, Cohen and Levin |101| by arguments based on the Heisenberg uncer- 
tainty relation, and independently with the help of virial theorem appl ied for a QCD 
motivated confining potential. See Fig. p|, r eproduced from Ref. |102]. Note that a 
similar decrease was predicted in Ref. [103|, which would depend not on the mass, 
but on the transverse mass of the particles, if the particle production happens so 
that the position of the emission is very strongly correlated with the momentum 
of the emitted particle |103]. So, it would be timely to check whether the effect 
depends on the particle mass, or on the transverse mass. Although the side ra- 
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dius components indicate such a decrease in case of pions, similar measurements for 
kaons and A-s would be indispensable to clarify the origin of the observed behavior. 

The question arises: can the effects be explained in a unified framework, 

that characterizes the hadronization process in e+e^ annihilation? An explanation 
of the rather small effective size of the source of the A-s seems to be a challenge for 
the Lund string model. 

The three-dimensional analysis of the NA22 data on h -f p reactions indicated 
a strong decrease of all the characteristic radii with increasing values of transverse 
momenta of the pair in the NA22 experiment |^. A decrease of the effective source 
sizes with increasing values of the transverse mass for a given kind of particle is 
seen in heavy ion collisions, similarly to effect i) in particle physics. The property 
iii), the decrease of the effective source size with the increase of the mass of the 
particle is seen in heavy ion physics and is explained in terms of hydrodynamical 
expansion, similarly to the explanation of effect i), see Figs |l^ and |l6| in Section [ij. 
Can one give a unified explanation of these similarities between results of particle 
interferomctry in e++ e^, h + p and heavy ion physics? We do not yet know the 
answer to this question. 



7. Invariant Buda— Lund Particle Interferometry 

The n-particle Bose-Einstein correlation function of Eq. (|^) is defined as the ratio 
of the 71-particle invariant momentum momentum distribution divided by an n- 
fold product of the single-particle invariant momentum distributions. Hence these 
correlation functions are boost-invariant. 

The invariant Buda-Lund parameterization (or BL in short) deals with a boost- 
invariant, multi-dimensional characterization of the building blocks {ata-k^) of ar- 
bitrary high order Bose-Einstein correlation functions, based on Eqs (p|p^). The BL 
parameterization was developed by the Buda-Lund Collaboration in Refs [jisl,^. 

The essential part of the BL is an invariant decomposition of the relative mo- 
mentum q in the exp(ig • Aa;) factor into a temporal, a longitudinal and two trans- 
verse relative momentum components. This decomposition is obtained with the 
help of a time-like vector in the coordinate space, that characterizes the center of 
particle emission in space-time, see Fig. 

Although the BL parameterization was introduced in Ref. ||l8| for high energy 
heavy ion reactions, it can be used for other physical situations as well, where 
a dominant direction of an approximate boost-invariant expansion of the particle 
emitting source can be identified and taken as the longitudinal direction Tz- For 
example, such a direction is the thrust axis of single jets or of back-to-back two-jet 
events in case of high energy particle physics. For longitudinally almost boost- 
invariant systems, it is advantageous to introduce the boost invariant variable r 
and the space-time rapidity rj, 



r/ ^ 0.51og[(t + r,)/(t-r,)] 



(79) 
(80) 



28 



T. Csorgo 




Z 



Fig. 10. Space-time picture of particle emission for a given fixed mean momentum 
of the pair. The mean value of the proper-time and the space-time rapidity distri- 
butions is denoted by r and rj. As the rapidity of the produced particles changes 
from the target rapidity to the projectile rapidity the \T(jj),rj{y)] variables scan the 
surface of mean particle production in the (t, r^) plane. 



Similarly, in momentum space one introduces the transverse mass mt and the ra- 
pidity y as 



mt = VE^-pI (81) 
y = 0.5\og[{E + p,)/{E-p,)]. (82) 

The source of particles is characterized in the boost invariant variables r, mt and rj— 
y. For systems that are only approximately boost-invariant, the emission function 
may also depend on the deviation from mid-rapidity, j/o- The scale on which the 
approximate boost-invariance breaks down is denoted by A77, a parameter that is 
related to the width of the rapidity distribution. 

The correlation function is defined with the help of the Wigner function formal- 
ism, Eq. (p^), the intercept parameter A* is introduced in the core/halo picture of 
Eq. (^) . The case of n = 2 particles and a chaotic core with Pc = was discussed 
in Ref. [|l8|. In the following, we evaluate the building block for arbitrary high 
order Bose-Einstein correlation functions. We assume for simplicity that the core 
is fully incoherent, Pc{j) = in Eq. (|6^). A further simplification is obtained if we 
assume that the emission function of Eqs factorizes as a product of an ef- 

fective proper-time distribution, a space-time rapidity distribution and a transverse 
coordinate distribution |104, 1^: 



Sc{x,K)d'^x — H^{T)Gr{r])I^{rx,ry) dr Tdr]drxdry. (83) 



The subscript * stands for a dependence on the mean momentum K, the mid- 
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rapidity yo and the scale of violation of boost-invariance A77, using the symbolic 
notation = f[K, yo, A77]. The function H^(t) stands for such an effective proper- 
time distribution (that includes, by definition, an extra factor r from the Jacobian 
d*x — dT T drj^drxdry, in order to relate the two-particle Bose-Einstein correlation 
function to a Fourier transformation of a distribution function in r). The effective 
space-time rapidity distribution is denoted by G, (77), while the effective transverse 
distribution is denoted by /*(rj^,rj^). In Eq. (|83|), the mean value of the proper-time 
r is factored out, to keep the distribution functions dimensionless. Such a pattern 
of particle production is visualized in Fig. |lO|. 

In case of hydrodynamical models, as well as in case of a decaying Lund 



strings |104,20|, production of particles with a given momentum rapidity y is lim- 
ited to a narrow region in space-time around rj and r. If the sizes of the effective 
source are sufficiently small (if the Bose-Einstein correlation function is sufficiently 
broad), the exp(i(7 • Ax) factor of the Fourier transformation is decomposed in the 
shaded region in Fig. as 

exp[i(g°Ai- g^Ar^)] ~ exp[i(Q=AT - QhtAt;)], (84) 
exp[-i((jj;Arj; + igj^Ary)] = exp[-i((3: Ar: + Q.. Ar..)]. (85) 

The invariant temporal^ parallel^ sideward, owtward (and perpendicular) relative 
momentum components are defined, respectively, as 

Q= = gocosh[^] - g^sinh[77], (86) 
Qii = cosh[7j] - 90 sinh[77], (87) 



0.. - {q,Ky-qyK,)/^Kl + Kl (88) 
Q: = {q,K,+qyKy)/^Kl + Kl (89) 
Qi- = Jql + ql = ^Q^ + Q^, (90) 



The time-like normal- vector n indicates an invariant direction of the source in co- 
ordinate space ||l^. It is parameterized as = (cosh[77], 0, 0, sinh[77]), where 7/ is a 
mean space-time rapidity ||l^,^,|2^ . The parameter is one of the fitted parameters 
in the BL type of decomposition of the relative momenta. The above equations are 
invariant, they can be evaluated in any frame. To simplify the presentation, in the 
following we evaluate q and rj in the LCMS. The acronym LCMS stands for the Lon- 
gitudinal Center of Mass System, where the mean momentum of a particle pair has 
vanishing longitudinal component, = 0.5(fci_2 + k2,z) = 0. In this frame, intro- 



duced in Ref. |f04|, K is orthogonal to the beam axis, and the time- like information 
on the duration of the particle emission couples to the out direction. The rapidity 
of the LCMS frame can be easily found from the measurement of the momentum 
vectors of the particles. As rj is from now on a space-time rapidity measured in the 
LCMS frame, it is invariant to longitudinal boosts: rj' = {rj — y) — {0 — y) = rj. 

The symbolic notation for the side direction is two dots side by side as in Q... 
The remaining transverse direction, the out direction was indexed as in Q,, in an 



30 



T. Csorgo 



attempt to help to distinguish the zeroth component of the relative momentum Qq 
from the out component of the relative momentum Q, = Qo ~ Qout, Qo 7^ Qo- 
Hence K- — |K_l| and K = 0. The geometrical idea behind this notation is 
explained in details in Ref. |^ . 

The perpendicular (or transverse) component of the relative momentum is de- 
noted by Q±. By definition, Q.., Q- and Q± are invariants to longitudinal boosts, 
and = -q-q = Q^+Q^ + Qf, - Ql. 

With the help of the small source size (or large relative momentum) expansion 
of Eq. (p^, the amplitude Sc( 1,2) = s*(l,2) that determines the arbitrary order 
Bose-Einstein correlation functions in Eq. (^3|) can be written as follows: 

,^ 2) = H4Q^)G4Q^^)l{Q.., Q..) 

^,(0)G40)/,(G,0) ■ ^ ' 

This expression and Eq. (|63| ) yield a general, invariant, multi-dimensional Buda- 
Lund parameterization of order n Bose-Einstein correlation functions, valid for all 
n. The Fourier-transformed distributions are defined as 



i?,(Q=) = / drexp(iQ=r)i/4r), (92) 
Jq 

/oo 
d^?exp(-^Q||r?7)G,(7?), (93) 
-00 

/OO pOO 
dr., dr„exp{-iQ.r.,-iQ,r„)h{r.,,r„). (94) 
~oo J — 00 

As a particular case of Eqs ( |9l| , |63| ) for n = 2 and Pc{j) = 0, the two-particle 
BECF can be written into a factorized Buda-Lund form as 

iHADt iG.m I'.Mf 
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Thus, the BL results are rather generic. For example, BL parameterization may in 
particular limiting cases yield the power-law, the exponential, the double-Gaussian, 
the Gaussian, or the less familiar oscillating forms of Eq. (128), see also Ref. [ p7| . 
The Edgeworth, the Laguerre or other similarly constructed low-momentum ex- 
pansions can be applied to any of the factors of one variable in Eq. ( p5| ) to 
characterize these unknown shapes in a really model-independent manner, relying 
only on the convergence properties of expansions in terms of complete orthonormal 
sets of functions [ |65| . 

In a Gaussian approximation and assuming that R, — R.. = the Buda- 
Lund form of the Bose-Einstein correlation function reads as follows: 



C2(ki, ka) =1+A, cxp i-RlQl - R\Q\ - RIQ 



(96) 



where the 5 fit parameters are A*, R=, R\\, R± and the value of ry that enters the 
definitions of Q= ^^^d Q\\ in Eqs ( |86| , |87| ). The fit parameter R= reads as _R-time- 
like, and this variable measures a width of the proper-time distribution H^,. The 
fit parameter _R|| reads as i?-parallel, it measures an invariant length parallel to 
the direction of the expansion. The fit parameter R± reads as i?-perpendicular or 
i?-perp. For cylindrically symmetric sources, R± measures a transversal rms radius 
of the particle emitting source. 



The BL radius parameters characterize the lengths of homogeneity |105l in a 
longitudinally boost-invariant manner. The lengths of homogeneity are generally 
smaller than the momentum-integrated, total extension of the source, they measure 
a region in space and time, where particle pairs with a given mean momentum K 
are emitted from. 

The following Edgeworth expansion can be utilized to characterize non-Gaussian 
multidimensional Bose-Einstein correlation functions, in a longitudinally boost- 
invariant manner: 



C2(ki,k2) 



1 + Xe exp{-QiRl 



l|..||-Qii?i)x 
1 + ^H3{V2Q=R=) + !l^Hi{V2Q=R=) - 
1 + ^i?3(\/2Q||i?||) + ^i/4(V2Q||i?||) + 
l + ^H,{V2Q^R^) + ^ 



Hi{V2Q^RA 



(97) 



This yields 5 free scale parameters for cylindrically symmetric, longitudinally ex- 
panding sources, and three series of shape parameters. The scale parameters are 
A^;, R=, R± and rj, that characterize the effective source at a given mean 
momentum, by giving the vertical scale of the correlations, the invariant tempo- 
ral, longitudinal and transverse extensions of the source and its invariant direction, 
which is the space-time rapi dity of the effective source in the LCMS frame (the 
frame where ki^^ + ^2.^ — 0, [104|). The three series of shape parameters are ^3=, 
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K4=, ... , K311 , K411 , ... , K3±, K4:A_, ... . Each of these parameters may depend on the 
mean momentum K. 

A multi-dimensional Laguerre or a mixed Edgeworth-Laguerre expansion can 
be introduced in a similar manner by replacing the Edgeworth expansion in Eq. ( |9^ ) 
by a Laguerre one in any of the principal directions. 

In Eqs ( |9^ , p6| ), the spatial information about the source distribution in (rx^Vy) 
was combined to a single perp radius parameter R±. In a more general Gaussian 
form, suitable for studying rings of fire and opacity effects, the Buda-Lund invariant 
BECF can be denoted as 



C2(ki, ka) - 1 + A* exp i^-RtQi - i?jfQjj - R^Qt - RfQ.) ■ (98) 

The 6 fit parameters are A», R=, i?.., R-, and 77, all are in principle functions of 
{'KjHq, Arj). Note that this equation is identical to Eq. (44) of Ref. iQ, rewritten 
into the new, symbolic notation of the Lorentz-invariant directional decomposition. 

The above equation may be relevant for a study of expanding shells, or rings of 
fire, as discussed first in Ref. We shall argue, based on a simultaneous analysis 
of particle spectra and correlations, and on recently found exact solutions of non- 
relativistic fireball hydrodynamics that an expanding, spherical shell of fire is 
formed protons in 30 AMeV ''"Ar -I^^^Au reactions, and that a two-dimensional, 
expanding ring of fire is formed in the transverse plane in NA22 h + p reactions at 
CERN SPS. The experimental signatures for the formation of these patterns will 



be discussed in Section 11 



Opacity effects, as suggested recently by H. Heiselberg |ll06| , also require the 
distinction between R,, and R.. The lack of transparency in the source may result 
in an effective source function, that looks like a crescent in the side-out reference 



frame [106|. When integrated over the direction of the mean momentum, the effec- 
tive source looks like a ring of fire in the {r^.Ty) frame. 

The price of the invariant decomposition of the basic building blocks of any 
order Bose-Einstein correlation functions in the BL parameterization is that the 
correlation functions cannot be directly binned in the BL variables, as these can 
determined after the parameter 77 is fitted to the data — so the correlation function 
has to be binned first in some directly measurable relative momentum components, 
e.g. the (side, out, long) relative momenta in the LCMS frame, as discussed in the 
next subsection. After fitting rj in an arbitrary frame, the BECF can be re-binned 
into the BL form. 



7. 1 . Gaussian parameterizations of BE Correlations 

We briefly summarize here the Bertsch-Pratt and the Yano-Koonin parameteriza- 
tion of the Bose-Einstein correlation functions, to point out some of their advantages 
as well as drawbacks and to form a basis for comparison. 
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7.1.1. The Bertsch-Pratt parameterization 

The Bertsch-Pratt (BP) parameterization of Bose-Einstein correlation functions is 
one of the oldest, widely used multi-dimensional decomposition, called also as the 
side-out -longitudinal decomposition p8|,p7|. 

This directional decomposition was devised to extract the contribution of a 
long duration of particle emission from an evaporating Quark-Gluon Plasma, as 
expected in the mixture of a hadronic and a QGP phase if the re-hadronization 
phase transition is a strong first order transition. 

The BP parameterization in a compact form reads as 

C2(ki,k2) = l + Xcxp[-RlQl-RlQl-R^Qf-2RliQiQo]. (99) 

Here index o stands for out (and not the temporal direction), s for side and I 
for longitudinal. The out-longitudinal cross-term was introduced by Chapman, 
Scotto and Heinz in Refs |^,^ — this term is non- vanishing for axially symmetric 
systems, if the source is not fully boost-invariant, or if the measurement is made 
not at mid-rapidity. In a more detailed form, the mean momentum dependence of 
the various components is shown as 

C2(ki,k2) = l + A(K)exp[-i?2(K)Q2(K)-i?2(K)g2(K) 

-i?2(K)Q2 _ 2Rli(K)QiQ,iK)] , (100) 

where the mean and the relative momenta are defined as 

K = 0.5(ki+k2), (101) 

Ak = ki-k2, (102) 

Qi = k,^i-k,a, (103) 

Qo = Qo(K) = Ak-K/|K|, (104) 

Qs - Qs{K) = |AkxK|/|K|. (105) 

It is emphasized that the BP radius parameters are also measuring lengths of ho- 



mogeneity 1 105]. Not only the radius parameters but also the decomposition of the 
relative momentum to the side and the out components depends on the (direction 
of) mean momentum K. 

In an arbitrary frame, Gaussian radius parameters can be defined, and some- 
times they are also referred to as BP radii, when the spatial components of the 
relative momentum vector are taken as independent variables. The BP radii reflect 
space-time variances |2^, ^ of the core ||7|] of the particle emission, if a Gaussian 
approximation to the core is warranted: 

C2(ki,k2) = 1 + A,(K) exp(-i?^(K)Ak,Ak,), (106) 
A4K) = [7Vc(K)/7V(K)]2, (107) 
RliK) = (£,£,)c-(£,)c(x,)c, (108) 
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= Xi-iiit, (109) 
(/(x,k))c = y" d^a;/(:r,k)5c(.T,k)/ ^ d'^xS^^x,'^), (110) 

where S^ix^ k) is the emission function that characterizes the central core and sub- 
scripts i or J stand for x, y or z, i.e. any of the spatial directions in the frame of 
the analysis. This method is frequently called as "model-independent" formulation, 
because the applied Gaussian approximation is independent of the functional form 
of the emission function S'(X,k) [ p^ . In the literature, this result is often over- 
stated, it is claimed that such a Taylor expansion would provide a general "proof 
that multi-dimensional Bose-Einstein correlation functions must be Gaussians. Al- 
though the "proof is indeed not depending on the exact shape of S{x,K), it rehes 
on a second order Taylor expansion of the shape of the correlation function around 
its exact value at Q = 0. At this point not only the derivatives of the correlation 
function are unmeasurable, but the very value of the correlation function C2(0) is 
unmeasurable as well, see Figs |4| and § for graphical illustration. For exponential 
or for power-law type correlations, the building block Sc{q,K) of the correlation 
function is not analytic at Q = 0, so a Taylor expansion cannot be applied in their 
case. For the oscillatory type of correlation functions, the Gaussian provides a good 
approximation in the experimentally unresolvable low Q domain, but it misses the 
structure of oscillations at large values of Q, which appear because S{x, K) has 
more than one maxima, like a source distribution of a binary star. Thus, the ex- 
act shapes of multi-dimensional BECF-s cannot he determined a priori and in case 
of non-Gaussian correlators one has to evaluate more (but still not fully) model- 
independent relationships, for example Eqs (|l^,^^l]), which are valid for broader 
than Gaussian classes of correlation functions. 

Note that the tails of the emission function are typically dominated by the halo 
of long-lived resonances Sh{x, k) and even a small admixture of e.g. r] and 77' mesons 
increases drastically the space-time variances of particle production, and makes the 
interpretation of the BP radii in terms of space-time variances of the total emission 
function S — Sc + Sh unreliable both qualitatively and quantitatively, as pointed 
out already in Ref. f?^ . 



In the Longitudinal Center of Mass System (LCMS, Ref. \104\), the BP radii 



have a particularly simple form |104], if the coupling between the r^ and the t 



coordinates is also negligible, (fxt) = {fx){t)- 

RUK) = (r~2),, (111) 

RliK) = (rt)c + A'(P)c, (112) 

i?2(K) = (r1)e, (113) 

RliiK) = -(3ti))c, (114) 

where x — x — (x) . Although this method cannot be applied to characterize non- 
Gaussian correlation functions, the the above form has a number of advantages: 
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it is straightforward to obtain and it is easy to implement for a numerical evaluation 
of the BP radii of Gaussian correlation functions . 

In the LCMS frame, information on the duration of the particle emission couples 
only to the out direction. This is one of the advantages of the LCMS frame. Using 
the BP, the time distribution enters the out radius component as well as the out- 
long cross-term. Other possible cross-terms were shown to vanish for cylindrically 
symmetric sources [p2|,p3|. 

For completeness, we give the relationship between the invariant BL radii and 
the BP radii measured in the LCMS, if the BL forms are given in the Gaussian 
approximation of Eq. (Bq) : 



i?2 = (115) 

i?2 ^ i?2 +/32[cosh2(77)i?2^ +sinh2(77)i?|], (116) 

Rli = -Ptsmh{r])cosh{r]){Ri+ Rj), (117) 

Rf cosh2(7y)i?| +sinh^(^)i?L (118) 

where the dependence of the fit parameters on the value of the mean momentum, 
K is suppressed. The advantage of the BP parameterization is that there are no 
kinematic constraints between the side, out and long components of the relative 
momenta, hence the BP radii are not too difficult to determine experimentally. A 
drawback is that the BP radii are not invariant, they depend on the frame where 
they are evaluated. The BP radii transform as a well-defined mixture of the invariant 
temporal, longitudinal and transverse BL radii, given e.g. in Ref. [T^ . 

7.1.2. The Yano-Koonin-Podgoretskii parameterization 

A covariant parameterization of two-particle correlations has been worked out for 



non-expanding sources by Yano, Koonin and Podgoretskii (YKP) |107, 108 1. This 
parameterization was recently applied to expanding sources by the Regensburg 
group | 109| , 110 1, by allowing the YKP radius and velocity parameters be momentum 



dependent: 

C2(ki,k2) = l + Cxp\-Rl{K)ql-Rj{K){ql-ql) 



Rl{K) + R^^{K)){q-U{K)f , (119) 



where the fit parameter J7(K) is interpreted 109, 110| as a four-velocity of a fluid- 



element I 111]. (Note that in YKP index refers to the time- like components). This 
generalized YKP parameterization was introduced to create a diagonal Gaussian 
form in the "rest frame of a fluid-element" . 

This form has an advantage as compared to the BP parameterization: the three 
extracted YKP radius parameters, R±, and Rq are invariant, independent of the 
frame where the analysis is performed, while transforms as a four-vector. The 
price one has to pay for this advantage is that the kinematic region may become 
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rather small in the go, qi, qi_ space, where the parameters are to be fitted, as follows 
from the inequalities = —q ■ q> and g§ > 0: 



o<qo< il 



(120) 



and the narrowing of the regions in q^ — q^ with decreasing q±^ makes the experi- 
mental determination of the YKP parameters difficult, especially when the analysis 
is performed far from the LCMS rapidities [or more precisely from the frame where 
f/'' = (1,0,0,0)]. 

Theoretical problems with the YKP parameterization are explained as follows. 
a) The YKP radii contain components proportional to l//3t, which lead to divergent 
terms for particles with very low pt |10£, L10[^ h) The YKP fit parameters are not 



even defined for all Gaussian sources |10£,110|. Especially, for opaque sources, for 



expanding shells, or for rings of fire with (f^) < (r^) the algebraic relations defining 
the YKP "velocity" parameter become ill-defined and result in imaginary values 
of the YKP "velocity", |lO§,[llo|. c) The YKP "flow velocity" C/^(K) is defined 
in terms of space-time variances at fixed mean momentum of the particle pairs 
[ 109 , 110 [, corresponding to a weighted average of particle coordinates. In contrast, 
the local flow velocity u^{x) is defined as a local average of particle momenta. 
Hence, in general [/^(K) ^ u^^{x), and the interpretation of the YKP parameter 
t/^(K) as a local flow velocity of a fluid does not correspond to the principles of 
kinetic theory. 



8. Hydrodynamical Parameterization a la Buda— Lund (BL-H) 

The Buda-Lund hydro parameterization (BL-H) was invented in the same paper 
as the BL parameterization of the Bose-Einstein correlation functions jl^, but 
in principle the general BL forms of the correlation function do not depend on 
the hydrodynamical ansatz (BL-H). The BL form of the correlation function can 
be evaluated for any, non-thermalized expanding sources, e.g. for the Lund string 
model also. 

The BL-H assumes, that the core emission function is characterized with a 
locally thermalized, volume-emitting source: 

S,{x,k)d^x = fc^rf^S^l^) ^ (121) 

The degeneracy factor is denoted by g, the four- velocity field is denoted by u^{x), 
the temperature field is denoted by T{x), the chemical potential distribution by ii{x) 
and s = 0, —1 or 1 for Boltzmann, Bose-Einstein or Fermi-Dirac statistics. The 
particle flux over the freeze-out layers is given by a generalized Cooper-Frye factor, 
assuming that the freeze-out hypersurface depends parametrically on the freeze-out 
time r and that the probability to freeze-out at a certain value is proportional to 
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H{r), 



k^d'^T,^{x) = rrit cosh[ri — y]H{T)dT rdr] dry. 



(122) 



The four- velocity u^(x) of the expanding matter is assumed to be a scaling longi- 
tudinal Bjorken flow appended with a linear transverse flow, characterized by its 
mean value (ut), see Refs |lS,23, 29|: 



{x) = i cosh[77] cosh[77t], sinh[?7t] — , sinh[r]t] — , sinh[r]] cosh[7yt] 



n 



n 



sinh[r/f 



{ut)rt/RG, 



(123) 



r^]^/^. Such a flow profile, with a time-dependent radius parameter 



with rt = [r. 

Rg, was recently shown to be an exact solution of the equations of relativistic 
hydrodynamics of a perfect fluid at a vanishing speed of sound, Ref. |40j . 

Instead of applying an exact hydrodynamical solution with evaporation terms, 
the BL-H characterizes the local temperature, flow and chemical potential distribu- 
tions of a cylindrically symmetric, finite hydrodynamically expanding system with 
the means and the variances of these distributions. The hydrodynamical variables 
1/T(x), iJ,{x)/T{x), are parameterized as 



^i{x) 


Mo 


2 , 2 


T{x) 


To 


2Rl 


1 


1 




T{x) 






AT 



\ T /t 2A7 



(124) 
(125) 



the temporal distribution of particle evaporation H{t) is assumed to have the form 
of 



2A7 



(126) 



(27rAT2)3/2 

and it is assumed that the widths of the particle emitting sources, e.g. Rg and Arj do 
not change significantly during the course of the emission of the observable particles. 
The parameters {AT/T)r and {AT/T)t control the transversal and the temporal 
changes of the local temperature profile, see Refs |2^,|l9|,[lJ] for further details. This 
formulation of the BL hydro source includes a competition between the transversal 
flow and the transverse temperature gradient, in an analytically tractable form. 
In the analytic evaluation of this model, it is assumed that the transverse flow is 
non-relativistic at the point of maximum emissivity [23| , the temperature gradients 
were introduced following the suggestion of Akkelin and Sinyukov |112|. 

Note that the shape of the profile function in -q is assumed to be a Gaussian in 
Eq. (124) in the spirit of introducing only means and variances. However, in Ref. fl^ 
a formula was given, that allows the reconstruction of this part of the emission 
function from the measured double-differential invariant momentum distribution in 
a general manner, for arbitrary sources with scaling longitudinal expansions. 
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8.1. Correlations and spectra for the BL-Hydro 

Using the binary source formulation, reviewed in the next section, the invariant 
single particle spectrum is obtained as 

Ni{k) = —^EVC 7 r . (127) 



exp 



Tix) T{x) 



The two-particle Bose-Einstein correlation function was evaluated in the binary 
source formalism in Ref. p7[ |: 

C2(ki,k2) = l + A,r!(Q||) exp (-Q|i?5 - Qii?L - Qii?l) , (128) 

where the pre-factor f2((3||) induces oscillations within the Gaussian envelope as a 
function of Q\\. This oscillating pre-factor satisfies < ^{Q\\) < 1 and f7(0) = 1. 
This factor is given as 

n{Q\\) = cos2(Q[|i?|| A77) + sin2(g||i?[| A77)tanh^(77). (129) 

The invariant BL decomposition of the relative momentum is utilized to present the 
correlation function in the simplest possible form. Although the shape of the BECF 
is non-Gaussian, because the factor ^{Q\\) results in oscillations of the correlator, 
the result is still explicitly boost-invariant. Although the source is assumed to 
be cylindrically symmetric, we have 6 free fit parameters in this BL form of the 
correlation function: A*, i?=, i?_L, rj and A?/. The latter controls the period of 
the oscillations in the correlation function, which in turn carries information on the 
separation of the effective binary sources. This emphasizes the importance of the 
oscillating factor in the BL Bosc-Einstcin correlation function. 

The parameters of the spectrum and the correlation function are the same, de- 
fined as follows. In the above equations, a means a momentum-dependent average 
of the quantity a. The average value of the space-time four-vector x is parame- 
terized by (T^r],r^,T\j), denoting longitudinal proper-time, space-time rapidity and 
transverse directions. These values are obtained in terms of the BL-H parameters 
in a linearized solution of the saddle-point equations as 

T = To, (130) 

V = {yo-y)/[l + ^V^nit/To], (131) 

= (ut)RG = — , (132) 

= 0. (133) 



In Eq. ( 127 ) , E stands for an average energy, Y for an average volume of the effective 



source of particles with a given momentum k and C for a correction factor, each 
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defined in the LCMS frame: 

E = mt cosh(77), (134) 

V = (27r)ii?|| i?l (135) 
At 

C = exp(A^V2) /v/a7. (136) 

The average invariant volume V is given as a time-averaged product of the transverse 
area R± and the invariant longitudinal source size i?|| , given as 





= Rl = R^ ^ Rll [I + {{utf + {AT/T)r) E/To] , 


(137) 




= r^Af, 


(138) 




= Ariy {1 + Afj^E/To), 


(139) 


rI 


= At^ = ArV (l + {AT/T)tE/TQ). 


(140) 



This completes the specification of the shape of particle spectrum and that of the 
two-particle Bose-Einstein correlation fmiction. These results for the spectrum 
correspond to the equations given in Ref. jl^] although they are expressed here 
using an improved notation. 

In a generali zed form, the thermal scales are defined as the E/Tq — )■ oo limit 
of Eqs (137- 140| ), while the geometrical scales correspond to dominant terms in 



the E/To —> limit of these equations. In all directions, including the temporal 
one, the length-scales measured by the Bose-Einstcin correlation function arc dom- 
inated by the smaller of the thermal and the geometrical length-scales. As shown 



in Sections 11 and |1^, the width of the rapidity distribution and the slope of the 
transverse-mass distribution is dominated by the bigger of the geometrical and the 
thermal length-scales. This is the analytic reason, why the geometrical source sizes, 
the flow and temperature profiles of the source can only be reconstructed with 
the help of a simultaneous analysis of the two-particle Bose-Einstein correlation 
functions and the single-particle momentum distribution JT6|-pO|. 

If the geometrical contributions to the HBT radii arc sufficiently large as com- 
pared to the thermal scales, they cancel from the measured HBT radius parameters. 
In this case, even if the geometrical source distribution for different particles (pions, 
kaons, protons) were different, the HBT radii (lengths of homogeneity) approach 
a scaling function in the large E/Tq limit. Up to the leading order calculation 
in the transverse coordinate of the saddle-point, this model predicts a scaling in 
terms of E, which variable coincides with the transverse mass mt at mid-rapidity. 
Phenomenologically, the scaling law can be summarized as Ri oc m"' , where i in- 
dexes the directional dependence, and the exponent may be slightly rapidity 
dependent, due to the difference between E and mt, and it may phenomenologi- 
cally reflect the effects of finite size corrections as well. Note also that such a scaling 
limiting case is only a possibility in the BL-H, valid in certain domain of parameter 



40 



T. Csorgo 



space, but it is not a necessity. The analysis of Pb + Pb collisions at 158 AGeV 
indicates that BL-H describes the data fairly well, but the longitudinal radius com- 
ponent exhibits different scaling behavior from the transverse radii, see Section ^ 
for more details. 



9. Binary Source Formalism 

Let us first consider the binary source representation of the BL-H model. The 
two-particle Bose-Einstein correlation function was evaluated in Ref. ||l8| only in a 
Gaussian approximation, without applying the binary source formulation. An im- 
proved calculation was recently presented in Ref. , where the correlation function 
was evaluated using in the binary source formulation, and the corresponding oscil- 
lations were found. 

Using the exponential form of the cosh[?7 — y] factor, the BL-H emission function 
Sc{x, k) can be written as a sum of two terms: 



Sc{x,k) = 0.5[S'+(x,k)-HS'_(x,k)], (141) 
-mtexp[±7y=F y]j 



S±{x,k) = -f-^mtexp[±7^Ty]H4T)————^, (142) 
(27r)J [fB[x,k)+s] 



/B(a;,k) = exp 



T{x) 



(143) 



Let us call this splitting as the binary source formulation of the BL-H parame- 
terization. The effective emission function components are both subject to Fourier 
transformation in the BL approach. In an improved saddle-point approximation, the 
two components S+(x,k) and S-{x,k) can be Fourier-transformed independently, 
finding the separate maxima (saddle point) x+ and X- of S+{x,k) and S-{x,k), 
and performing the analytic calculation for the two components separately. 

The oscillations in the correlation function are due to this effective separation of 
the pion source to two components, a splitting caused by the Cooper-Frye flux term. 
These oscillations in the intensity correlation function are similar to the oscillations 
in the intensity correlations of photons from binary stars in stellar astronomy |113[ . 

Due to the analytically found oscillations, the presented form of the BECF goes 
beyond the single Gaussian version of the saddle-point calculations of Refs [p2|,p3[. 
This result goes also beyond the results obtainable in the YKP or the BP param- 
eterizations. In principle, the binary-source saddle-point calculation gives more 
accurate analytic results than the numerical evaluation of space-time variances, as 
the binary-source calculation keeps non-Gaussian information on the detailed shape 
of the Bose-Einstein correlation function. 

Note that the oscillations are expected to be s mall in the BL-H picture, and 
the Gaussian remains a good approximation to Eq. ( |l28|) , but with modified radius 
parameters. 
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9. 1 . The general binary source formalism 

In the previous subsection, we have seen liow effective binary sources appear in the 
BL-H model in high energy physics. However, binary sources appear generally: in 
astrophysics, in form of binary stars, in particle physics, in form of W"'"W~ pairs, 
that separate before they decay to hadrons. 

Let us consider first the simplest possible example, to see how the binary sources 
result in oscillations in the Bose-Einstein or Fermi-Dirac correlation function. Sup- 
pose a source distribution s{x — x+) describes e.g. a Gaussian source, centered on 
x+. Consider a binary system, where the emission happens from s-|_ = s{x — x-^-) 
with fraction /+, or from a displaced source, s_ — s{x — a;_), centered on x-, with 
a fraction /_ . For such a binary source, the amplitude of the emission is 

p{x) = f+s{x - x+) + f-s{x - X-), (144) 

and the normalization requires 

/+ + /- = !• (145) 
The two-particle Bose-Einstein or Fermi-Dirac correlation function is 

C{q) = l±|p(g)|2 = l±^g)|S(g)|2, (146) 

where + is for bosons, and — for fermions. The oscillating pre-factor Cl{q) satisfies 
< ^l{q) < I and r2(0) = 1. This factor is given as 

m = [{fl + f_)+2f+f_cos[qix+-x.)]]. (147) 

The strength of the oscillations is controlled by the relative strength of emission 
from the displaced sources and the period of the oscillations can be used to learn 
about the distance of the emitters. In the limit of one emitter (/+ = 1 and /_ = 0, 
or vice versa), the oscillations disappear. 

The oscillating part of the correlation function in high energy physics is expected 
to be much smaller, than that of binary stars in stellar astronomy. In particle 
physics, the effective separation between the sources can be estimated from the 
uncertainty relation to be x± = \x+ — a;_| w 2Ti/M^ k, 0.005 fm. Although this is 
much smaller, the effective size of the pion source, 1 fm, one has to keep in mind 
that the back-to-back momenta of the W+W^ pairs can be large, as compared to 
the pion mass. Due to this boost, pions with similar momentum may be emitted 
from different W-s with a separation which is already comparable to the 1 fm 
hadronization scale, and the resulting oscillations may become observable. 

In stellar astronomy, the separation between the binary stars is typically much 
larger than the diameter of the stars, hence the oscillations are well measurable. In 
principle, similar oscillations may provide a tool to measure the separation of the 
W+ from W~ in 4-jet events at LEP2. The scale of separation of W+W" pairs is a 
key observable to estimate in a quantum-mechanically correct manner the influence 
of the Bose-Einstein correlations on the reconstruction of the W mass. 
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In heavy ion physics, osciUations are seen in the long-range part of the p + p 



Fermi-Dirac correlation function [114|, with a half-period of Qh = 30 MeV. This 



implies a separation of a;± — nh/Qh ~ 20 fm, which can be attributed to interference 



between the the two peaks of the NA49 proton dn/dy distribution [115|, separated 
by Ay = 2.5. As for the protons we have Tq = 140 MeV, we can identify this 
rapidity difference with the space-time rapidity difference between the two peaks 
of the rapidity distribution. The longitudinal scale of the separation is then given 
by x± = 2t sinh(A7/p/2), which can be used to estimate the mean freeze-out time 
of protons, r = 7r?i/[2(5/i sinh(A77p/2)] w 6.4 fm/c, in a good agreement with the 
average value of r = 5.9 ± 0.6 as extracted from the simultaneous analysis of the 
single-particle spectra and HBT radii in NA44, NA49 and WA98 experiments in 
the Buda-Lund picture, as summarized in Section 

10. Particle Correlations and Spectra at 30 - 160 AMeV 

There are important qualitative differences between relativistic heavy ion collisions 
at CERN SPS and those at non-relativistic energies from the point of view of par- 
ticle sources. Low and intermediate energy reactions may create a very long-lived, 
evaporative source, with characteristic lifetimes of a few 100 fm/c, in contrast to 
the relatively short-lived systems of lifetimes of the order of 10 fm/c at CERN SPS. 
During such long evaporation times, cooling of the source is unavoidable and has to 
be included into the model. Furthermore, in the non-relativistic heavy ion collisions 
mostly protons and neutrons are emitted and they have much stronger final-state 
interactions than the pions dominating the final state at ultra-relativistic energies. 



see Refs [116-118| for recent reviews. 

The evolution of the particle emission in a heavy-ion collision at intermediate 
energies may roughly be described as: production of pre-equilibrium particles; ex- 
pansion and possible freeze-out of a compound source; possible evaporation from an 
excited residue of the source. Note though that this separation is not very distinct 
and there is an overlap between the different stages. The importance of the various 
stages above also depends on the beam energy and the impact parameter of the 



collision. See the review paper of Ref. |11S] for greater details 



Sophisticated microscopical transport descriptions ]119| ], such as the BUU 
(Boltzmann-Uehling-Uhlenbeck) and the QMD (Quantum Molecular Dynamics) 
models are well-known and believed to provide a reasonable picture of proton emis- 
sion in central heavy ion collisions from a few tenths up to hundreds of MeV per 
nucleon. However, the BUU model predicts too large correlations and underpredicts 
the number of protons emitted with low energies, for the reaction '^^Ar -I- ''^Sc at 



E = 120 and 160 MeV/nucleon, see Ref. [120|. This indicates that the simultane- 
ous description of two-particle correlations and single-particle spectra is a rather 
difficult task. For energies below a few tens of MeV per nucleon, where long-lived 
evaporative particle emission is expected to dominate, the measured two-proton 
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correlation functions were found to be consistent with compound-nucleus model 



predictions [121| 



A simultaneous analysis of proton and neutron single particle spectra and two- 
particle correlation was presented in Ref. This model calculation described 
the second stage above and, for long emission times, also part of the third stage. 
In Ref. p^ ], the competition among particle evaporation, temperature gradient 
and flow was investigated in a phenomenological manner, based on a simultaneous 
analysis of quantum statistical correlations and momentum distributions for a non- 
relativistic, spherically symmetric, three-dimensionally expanding, finite source. 
The model used can be considered as a non-relativistic, spherically symmetric ver- 
sion of the BL-H hydro parameterization j26j . 

The non-relativistic kinetic energy is denoted by Ek{h) = kV(2m). The fol- 
lowing result is obtained for the effective source size i?*: 

2 

" l + [{AT/T)rEk(k)+m{ut)^] /To ' ^^^^^ 

The analytic results for the momentum distribution and the quantum statistical 
correlation function are given in the Boltzmann approximation as 



Ni{k) = ^iJfe(k)K(k)exp 



(k-mu(r,(k)))2 /i(r,(k)) 



2mT(r,(k)) r(r,(k))_ 

(149) 

K(k) = [27ri?2(k)]'/' , (150) 
C{K,Ak) = l±cxp(-i?2(K)Ak^- At^Ai;^) . (151) 

The effects of final-state Coulomb and Yukawa interactions on the two-particle 
relative wave-functions are neglected in these analytic expressions. When comparing 
to data, the final-state interactions were taken into account, see Ref. for further 
details. 

These general results for the correlation function indicate structural similar- 
ity between the non-relativistic flows in low/intermediate energy heavy ion colli- 
sions ||l^,|2^, ^ and the transverse flow effects in relativistic high energy heavy 
ion and elementary particle induced reactions ||l9|, 20, Is!] . The radius parameters 



of the correlation function and the slopes of the single-particle spectra are momen- 
tum dependent both for the non-relativistic versions of the model, presented in 
Refs [|l^,|2^, ^ and for the model-class with scaling relativistic longitudinal flows, 
discussed in Refs ]l|-||,||l. 

Such a momentum-dependent effective source size has been seen in the proton- 
proton correlation functions in the ^^Al (^''N, pp) reactions at i? = 75 MeV/nucleon 



|117|: the larger the momentum of the protons the smaller the effective source size 
1 117 1, in qualitative agreement with Eq. ( |14^ ). 

This model was applied in Ref. to the reaction '*"Ar + ^^^Au at 30 MeV/ 
nucleon. With the parameter set presented in Table ^ we have obtained a simul- 
taneous description of the n and p single particle spectra as well as the nn and pp 
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Table 3. Parameter values obtained from fitting hydro parameters to n and p spec- 
tra and correlation functions, as measured by the CHIC Collaboration in 30 AMeV 
^"Ar + ^^^Au reactions 





Rg (fm) 


To (MeV) 


{AT/T)r 




Neutrons 


4.0 


3.0 


0.0 


0.018 


Protons 


4.0 


5.0 


0.16 


0.036 



correlation functions as given by Refs p2^~[l24[. See Ref. |26| for further details 
and discussions. 

The main effects of the temperature gradient are that it introduces i) a, momen- 
tum-dependent effective temperature which is decreasing for increasing momentum, 
resulting in a suppression at high momentum as compared to the Boltzmann dis- 
tribution; ii) a momentum-dependent effective source size which decreases with 
increasing total momentum. Agreement with the experimental data is obtained 
only if the time of duration of the particle emission was rather long, (t) w 520 fm 
with a variance of ~ 320 fm/c. 

The obtained parameter set reflects a moderately large system (Gaussian radius 
parameter Rq = 4.0 fm) at a moderate temperature {TQ{n) = 3 MeV and Tq{p) = 
5 MeV) and small flow. The neutrons and the protons seem to have different local 
temperature distributions: the neutron temperature distribution is homogeneous, 
while the temperature of the proton source decreases to Ts{p) = 4.3 MeV at the 
Gaussian radius, a difference that could be attributed to the difference between 
their Coulomb interactions An agreement between the model and the data 

was obtained only if some amount of flow was included | |2^ . 

After the completion of the data analysis, a new family of exact solutions of 
fireball hydrodynamics was found in Ref. , which features scaling radial Hubble 
flow, and an initial inhomogeneous, arbitrary temperature profile. The competition 
of the temperature gradients and flow effects were shown to lead to the formation 
of spherical shells of fire in this class of exact hydrodynamical solutions |^ , if the 
temperature gradient was stronger than the fiow, (AT/T)^ > m{ut)'^ /Tq. This is 
the case found from the analysis of proton spectra and correlations in Ref. [p6| , 
while the neutron data do not satisfy this condition. Assuming the validity of 
non-relativistic hydrodynamics to characterize this reaction, one finds that a slowly 
expanding, spherical shell of fire is formed by the protons, while the neutrons remain 
in a central, slightly colder and even slower expanding, normal fireball in 30 AMeV 
''"Ar + ""^^^Au heavy ion reactions. 
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11. Description of h + p Correlations and Spectra at CERN 
SPS 

The invariant spectra of tt^ mesons produced in (7r+/K+)p interactions at 250 
GeV/c are analysed in this section in the framework of the BL-H model of three- 
dimensionally expanding cylindrically symmetric finite systems, following the lines 
of Ref. The EHS/NA22 Collaboration has been the first to perform a detailed 
and combined analysis of single-particle spectra and two-particle Bose-Einstein cor- 
relations in high energy physics NA22 reported a detailed study of multi- 
dimensional Bose-Einstein correlations, by determining the side, out and the longi- 
tudinal radius components at two different values of the mean transversal momenta 
in (7r+/K+)p at CERN SPS energies |Q. It turned out, however, that the ex- 
perimental two-particle correlation data were equally well described by a static 
Kopylov-Podgoretskii parameterization as well as by the predictions of hydrody- 
namical parameterizations for longitudinally expanding, finite systems. In Refs 
we have shown, that the combined analysis of two-particle correlations and 
single-particle spectra may result in a dramatic enhancement of the selective power 
of data analysis. 



The double-differential invariant momentum distribution of Eq. (127) can be 
substantially simplified for one-dimensional slices ||l8 , 77 . 



i) At fixed mt, the rapidity distribution reduces to 

{y ~ yo? 



7Vi(k) = C™exp 



(152) 



2Ay2 

Ay2 = Aif+To/mt, (153) 

where Cm is an mt-dependent normalization coefficient and ?/o is defined above. The 
width parameter Ay^ extracted for difi^erent mj-slices is predicted to depend linearly 
on 1/mt, with slope Tq and intercept Ai]'^ . Observe, that this width is dominated 
by the bigger of the geometrical scale (A?/) and the thermal scale Tp/mt. 



Note that for static fireballs or spherically expanding shells (152) and (153) are 
satisfied with Arj = ||7^. Hence the experimental determination of the 1/mt de- 
pendence of the Ay parameter can be utilized to distinguish between longitudinally 
expanding finite systems versus static fireballs or spherically expanding shells. 

a) At fixed y, the to^ -distribution reduces to 

7Vi(k)-C,m^exp(^~^) , (154) 

where Cy is a y-dependent normalization coefficient and a is related to the effective 
dimensions of inhomogeneity in the source as a = 1 - deff/2 |l|]. The y-dependent 
"effective temperature" Tcs{y) reads as 

Tcff(y) = —-^ , (155) 

l + a(y-yo)^ 
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Fig. 11. The rapidity distributions of centrally produced pions (|y| < 1.5) for 
different mt-slices given. The curves are the fit results obtained analytically using 
the BL-H parameterization. 

where T* is the maximum of Tcsiy) achieved ai y — yo, and parameter a can be 
expressed with the help of the other fit parameters, see Refs Hl^. 

The slope parameter at mid-rapidity, is also determined by an interplay of 
the central temperature Tq the flow effects modeled by (u*)^ and the temperature 
difference between the surface and the center, as characterized by {AT/T)r |]l8| , p9| . 
Eq. (66) of Ref. can be rewritten as 



= To -4- m{uty 



To 



To + m(^)r 



(156) 
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Fig. 12. The m* distributions of centrally produced pions {\y\ < 1.5) for different 
y-slices given. The curves are the fit results obtained analytically using the BL-H 
parameterization. 



The approximations of Eqs ( 152 ) and (154) explicitly predict a specific narrowing of 
the rapidity and transverse-mass spectra with increasing and y, respectively (cf. 
(153) and (15£ )). The character of these variations is expected |0 to be different for 
the various scenarios of hadron matter evolution. These features of the spectra were 
found to be in agreement with the NA22 data |Q, and were utilized to reconstruct 
the particle source of h + p reactions in the {t,rz) plane. 
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Table 4. Fit results to NA22 h+p data at CERN SPS with a Buda-Lund hydro 
parameterization for \y\ < 1.5 



a Ary 


To (GeV) 


(ut) (AT/T), 


xVndf 


0.26 ±0.02 1.36 ±0.02 


0.140 ±0.003 


0.20 ±0.07 0.71 ±0.14 


642/683 



11.1. Combination with two-particle correlations 

As already mentioned in the introduction, more comprehensive information on geo- 
metrical and dynamical properties of the hadron matter evolution are expected from 
a combined consideration of two-particle correlations and single-particle inclusive 
spectra [||-|^,||,|ll|. 

At mid-rapidity, y = yo and in the LCMS where fci_z = — fc2,z, the effective BP 



radii can be approximately expressed form the BL-H parameterization as |18|: 



Rf = t^AtJ^ (157) 



Rl = rI+P^At^, (158) 
Rl = r]_ (159) 



with 



M, 



* (160) 



Arj^ A772 To 

= l + t((5+(¥^>.)' ^'''^ 

where parameters Atj^^Tq, (ut) and {AT/T)r are defined and estimated from the 
invariant spectra; Rq is related to the transverse geometrical rms radius of the 
source as i?G(rms) — V^Rg', t is the mean freeze-out (hadronization) time; At is 
related to the duration time At of pion emission and to the temporal inhomogeneity 
of the local temperature, as the relation At > At holds; the variable (3t is the 
transverse velocity of the pion pair. 

The effective longitudinal radius Ri, extracted for two different mass ranges, 
Mt = 0.26 ± 0.05 and 0.45 ± 0.09 GeV/c^ are found to be Ri = 0.93 ± 0.04 and 
0.70 ± 0.09 fm, resp ectively. This dependence on Mt matches well the predicted 
one. Using Eq. (|l59|) with Tq = 140±3 MeV and Ayy^ = 1.85 ±0.04 (Table |), one 



finds that the values of t extracted for the two different Mt regions are similar to 
each other: t = 1.44 ± 0.12 and 1.36 ± 0.23 fm/c. The averaged value of the mean 
freeze-out time is t = 1.4 ± 0.1 fm/c. 

The width of the (longitudinal) space-time rapidity distribution of the pion 
source was found to be Ai] = 1.36 ± 0.02. Since this value of A77 is significantly 
bigger than 0, the static fireballs or the spherically expanding shells fail to reproduce 
the NA22 single-particle spectra ||5^, although each of these models was able to 
describe the NA22 two-particle correlation data in Ref. [p9[. 
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The transverse-plane radii Ro and Rs were reported in Rcf. p3| for the whole 



Mt range are: Ro = 0.91 ± 0.08 fm and Rs = 0.54 ± 0.07 fm. Substituting in ( |l57| ) 
and dlS^ ), one obtains (at Pt = 0.484c [H): At = 1.3 ± 0.3 fm/c. The mean 
duration time of pion emission can be estimated as At > At = 1.3 ± 0.3 fm/c. 
A possible interpretation of At « t might be that the radiation process occurs 
during almost all the hydrodynamical evolution of the hadronic matter produced 
in meson-proton collisions. 



An estimation for the parameter Rq can be obtained from (158) and (16C) 
using the quoted values of Rg, Tq, (ut) and (AT/T) at the mean value of (Mt) = 

0. 31 ± 0.04 GeV/c (averaged over the whole Mt range): Rq = 0.88 ± 0.13 fm. The 
geometrical rms transverse radius of the hydrodynamical tube, i?G(™is) = \/2Rq = 

I. 2 ± 0.2 fm, turns out to be larger than the proton rms transverse radius. 

The data favour the pattern according to which the hadron matter under- 
goes predominantly longitudinal expansion and non-relativistic transverse expan- 
sion with mean transverse velocity (ut) — 0.20 ± 0.07, and is characterized by a 
large temperature inhomogeneity in the transverse direction: the extracted freeze- 
out temperature at the center of the tube and at the transverse rms radius are 
140 ± 3 MeV and 82 ± 7 MeV, respectively. 

II. 2. The space-time distribution of vr emission 

A reconstruction of the space-time distribution of pion emission points is shown in 
Fig. |l^, expressed as a function of the cms time variable t and the cms longitudinal 
coordinate z = r^. The momentum- integrated emission function along the z axis, 

1. e. at rt = {rx,ry) = (0,0) is given by 

Sit, z) oc exp r- exp f-^^) . (162) 



2At2 J \ 2Ary- 

It relates the parameters fitted to the NA22 single-particle spectrum and HBT 
radii to the particle production in space-time. The coordinates {t, z) are expressed 
with the help of the longitudinal proper-time t and space-time rapidity as 7] as 
(t cosh(77) , T sinh(?7) ) . 

We find a structure looking like a boomerang, i.e. particle production takes 
place close to the regions oi z — t and z = —t, with gradually decreasing probability 
for ever larger values of space-time rapidity. Although the mean proper-time for 
particle production is t = 1.4 fm/c, and the dispersion of particle production in 
space-time rapidity is rather small, Arj — 1.35 fm, we still see a characteristic long 
tail of particle emission on both sides of the light-cone, giving a total of 40 fm 
maximal longitudinal extension in z and a maximum of about 20 fm/c duration of 
particle production in the time variable t. 

In the transverse direction, only the rms width of the source can be directly 
inferred from the BP radii. However, the additional information from the analysis 
of the transverse momentum distribution on the values of {u)t and on the values of 
(AT/T)r can be used to reconstruct the details of the transverse density profile, as 
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Fig. 13. The reconstructed S{t, z) emission function in arbitrary units, as a func- 
tion of time t and longitudinal coordinate z. The best fit parameters of Arj = 1.36, 
yo = 0.082, At = 1.3 fm/c and t = 1.4 fm/c are used to obtain this plot. Note 
that before we made this reconstruction together with the NA22 Collaboration, 
only 1 fm^ area from this extended bumerang shape was visible to the intensity 
interferometry microscope. 

an exact, non-relativistic hydro solution was found in Rcf. | |39[ |, given in terms of 
the parameters {u)t and (AT/r)^ and using an ideal gas equation of state. Assum- 
ing the validity of this non-relativistic solution in the transverse direction, in the 
mid-rapidity range, one can reconstruct the detailed shape of the transverse density 
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NA22 (pi/K) + p at CERN SPS 




Fig. 14. The reconstructed S{rx,ry) emission function in arbitrary vertical units, 
as a function of the transverse coordinates Vx and ry. The shape has been recon- 
structed assuming the validity of a non-relativistic solution of hydrodynamics in the 
transverse direction, and using the values of Tq, (AT/T) and (ut) as obtained from 
the fits to the single-particle spectra. The root mean square width of the source 
distribution was obtained from the fits to the NA22 Bose-Einstein correlation func- 
tions. The momentum variables and the longitudinal and temporal variables are 
integrated over. 

profile. The result looks like a ring of fire in the {rx,ry) plane, see Fig. ^ In 
this hydro solution, {AT/T)r < m{u)'^/To corresponds to self-similar, expanding 
fireballs, while {AT/T)r > m{u)1/TQ corresponds to self-similar, expanding shells 
or rings of fire. 

Due to the strong surface cooling and the small amount of the transverse fiow, 
one finds that the particle emission in the transverse plane of h -I- p reactions at 
CERN SPS corresponds to a ring of fire. This transverse distribution, together with 
the scaling longitudinal expansion, creates an elongated, tube-like source in three 
dimensions, with the density of particle production being maximal on the surface 
of the tube. 
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12. Pb + Pb Correlations and Spectra at CERN SPS 



In Ref. |3^, an analysis similar to that of the NA22 Collaboration has been per- 
formed, with improved analytic approximations, using Fermi-Dirac or Bose-Einstein 
statistics (s = ±1) in the analytic expressions fitted to single particle spectra. The 
spectra were evaluated with the binary source method, the Bose-Einstein correla- 
tion functions were calculated with the saddle-point method without invoking the 
binary source picture. The analytical formulas for the BECF and IMD, as were used 
in the fits, were summarized in their presently most advanced form in Section ^, 
their development was described in Rcfs (l|-|o[ |3, |9|, . 
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Fig. 15. Result of simultaneous fits of the Buda-Lund hydro model to particle 
correlations and spectra in 158 AGeV Pb -I- Pb reactions at CERN SPS (data from 
the NA49 Collaboration) 

In case of homogeneous freeze-out temperatures, or particles with small masses. 
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Fig. 16. Result of simultaneous fits of the Buda-Lund hydro model to particle 
correlations and spectra in 158 AGeV Pb + Pb reactions at CERN SPS (data from 
the NA44 Collaboration) 



(163) 



Eq. (|156|) implies a linear rise of the slope with m [[18| as 



T,(m) = To + m{u)t, 



For heavy particles, or for large, non-vanishing temperature gradients, a flattening 
of the initial linear rise is obtained ll§l as 



r,(m) = To 



(164) 



This means that very heavy particles resolve the temperature inhomogeneities of 
the source, and they are produced with a mass-independent effective slope param- 
eter in the BL-H parameterization, if Ta/m becomes smaller than the temperature 
inhomogeneity. In a general case, the T^{'m) function starts with an initial linear 
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m dependence, with a slope given by the transverse flow {u)t, then T^,{m) flattenes 
out to a mass-independent value if the source has temperature inhomogeneities in 
the transverse direction. Such a behavior was reported by Pb + Pb heavy ion 
experiments at CERN SPS g. The central temperature is Tq w 140 MeV, 
the flattening of the slopes sets in at about m ~ 1400 MeV [p|7 which then leads 
to about 10% temperature inhomogeneity in the transverse direction of the Pb + 
Pb source. This estimate is in a good agreement with the results of the combined 
analysis of the single-particle spectra and the two-particle Bose-Einstein correlation 
functions, see Table 5. 



Table 5. Fit paramaters of Buda-Lund hydro (BL-H) in a simultaneous analysis 
of NA49, NA44 and preliminary WA98 spectra and correlation data 



Parameter 


NA49 


NA44 


WA98 


Averaged 


To [MeV] 


134 ±3 


145 ±3 


139 ±5 


139 ±6 


(ut) 


0.61 ±0.05 


0.57±0.12 


0.50 ±0.09 


0.55 ±0.06 


Rg [fm] 


7.3 ±0.3 


6.9± 1.1 


6.9 ±0.4 


7.1 ±0.2 


To [fm/c] 


6.1 ±0.2 


6.1 ±0.9 


5.2 ±0.3 


5.9 ±0.6 


At [fm/c] 


2.8 ±0.4 


0.01 ±2.2 


2.0 ± 1.9 


1.6 ± 1.5 


At] 


2.1 ±0.2 


2.4 ± 1.6 


1.7±0.1 


2.1 ±0.4 


{AT/T)r 


0.07 ±0.02 


0.08 ±0.08 


0.01 ±0.02 


0.06 ±0.05 


{AT/T)t 


0.16 ±0.05 


0.87 ±0.72 


0.74 ±0.08 


0.59 ±0.38 


xVndf 


163/98 = 1.66 


63/71 = 0.89 


115/108 = 1.06 


1.20 



The NA49, NA44 and WA98 data on single particle spectra of h~, identified tt, 
K and p as well as detailed rapidity and mt dependent HBT radius parameters are 
found to be consistent with each other as well as with BL-H. The BL-H fit results 
to these data sets is summarized in Table 5, Ref. [M. 



13. Comparison of h + p and Pb + Pb Final States at CERN 
SPS with Heavy Ion Reactions at Low and Intermediate 
Energies 

The final state of central Pb + Pb collisions at CERN SPS corresponds to a 
cylindrically symmetric, large {Rg = 7.1 ± 0.2 fm) and transversally homogenous 
(To = 139 ±6 MeV) fireball, expanding three-dimensionally with (ut) — 0.55 ±0.06. 
A large mean freeze-out time, t = 5.9 ± 0.6 is found with a relatively short du- 
ration of emission, At — 1.6 ± 1.5 fm, which is similar to the time-scale of emis- 
sion in the h ± p reaction. Note that the temporal cooling in Pb + Pb reac- 
tions seems to be stronger than in h ± p, which can be expained by the faster, 



Particle Interferomctry from 40 MeV to 40 TcV 



55 



three-dimensional expansion in the former case, as compared to the essentially one- 
dimensional expansion in the case of h -I- p reactions. By the time the particle 
production is over, the surface of Pb -I- Pb collisions cools down from 139 MeV to 
To/(l + (AT/To)^)/(l + {AT/To)t) ~ 83 MeV. It is very interesting to note that 
this value is similar to the surface temperature of Ts = 82 ± 7 MeV, found in h -|- 
p reactions as a consequence of the transverse temperature inhomogeneities, as de- 
scribed in Section ^ Ref. . Such snowballs with relatively low values of surface 
temperature Tg and a possible hotter core were reported first in 200 AGeV S -I- Pb 
reactions in Ref. ps| ]. 

Other hydro parameterizations, as reviewed in Ref. frequently neglect the 
effects of temperature inhomogeneities during the expansion and particle production 
stage. Energy conservation implies that the temperature cannot be exactly constant 
when particles are freezing out in a non-vanishing period of time from a three- 
dimensionally expanding source. 

The exact solution of non-relativistic, spherically symmetric fireball hydrody- 
namics implies [^8| , [l^ that Gaussian fireballs with spatially uniform temperature 
profiles satisfy the coUisionless Boltzmann equation. 

Fixing the temperature to a constant in the fits yields an average freeze-out 
temperature in the range of T/ = 110 ± 30 MeV §,|l|,|2|,|2|. 



Pb + Pb at CERN SPS 




Fig. 17. The reconstructed S{t, z) emission function in arbitrary units, as a func- 
tion of time t and longitudinal coordinate z, for 158 AGeV Pb -|- Pb reactions 
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Pb + Pb at CERN SPS 




S(x,y) 



Fig. 18. The reconstructed Sir^^ry) emission function in arbitrary units, as a 
function of the transverse coordinates and ry 

Based on the recently found new family of non-relativistic hydrodynamics 
and on the analysis of h + p single particle spectra and two-particle Bose-Einstein 
correlation function [Q, we concluded that the pion emission function S{rx,ry) 
in h + p reactions corresponds to the formation of a ring of fire in the transverse 
plane, because the transverse flow is rather small and because the sudden drop of 
the temperature in the transverse direction leads to large pressure gradients in the 
center and small pressure gradients and a density built-up at the expanding radius 
of the fire-ring. We presented arguments for a similar formation of a spherical shell 
of fire in the proton distributions at 30 AMeV ''"Ar -I- ^^^Au reactions. 

The formation of shells of fire seems thus to be of a rather generic nature, 
related to the initial conditions of self-similar radial flows. It is natural to ask the 
question: can we learn more about this phenomena in other physical systems? 

Radial expansion is a well established phenomena in heavy ion collisions from 



low energy to high energy reactions. See Refs |117, 118 for recent reviews and for 



example see Refs |125-12S] for the evidence of collective flow in central heavy ion 
collisions from 100 AMeV to 2 AGeV as measured by the FOPI Collaboration at 
GSI SIS. 

The FOPI Collaboration measured recently the proton-proton correlation func- 



tions at 1.93 AGeV Ni -|- Ni collisions |128[ . To interpret their data, they utilized a 
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Pb+PbatCERNSPS NA33(pi/KJ+patCERNSPS 




Fig. 19. Comparison of the reconstructed S{t,z) and S{rx,ry) emission functions 
for 250 GeV/c h + p reactions and for for 158 AGeV Pb + Pb reactions at CERN 
SPS. Note the different characteristic scales in the transverse and the temporal 
directions, and the different shapes of the transverse density distribution. 

version of the hydrodynamical solution, found in Ref. [Q. They assumed a linear 
flow profile, a Gaussian density distribution and a constant temperature. Such a so- 
lution of fireball hydrodynamics exists, but it corresponds to a collisionless Knudsen 
gas |38[^. A collisionless approximation has to break down. Indeed, only the peak 
of the FOPI proton-proton correlation function was reproduced by the collisionless 
model, however, the tails had to be excluded from the FOPI analysis. Perhaps it 
is worthwhile to search for a possible formation of shells of fire at the SIS energy 
domain, by re-analyzing the FOPI data 

14. Shells of Fire and Planetary Nebulae 

In transport calculations based on the Boltzmann-Uehling-Uhlenbeck equation, a 
formation of toroidal density distributions was predicted for central ^^Ar -|- '^^Sc 
collisions at i? = 80 AMeV in Ref. | ]120{ , which leads to ring-like configurations for 

S{rx, Ty). 

However, the clearest experimental observation of the development of expanding 
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Planetary nebula BD430: HST WFPC II Hp Image 



1 °| 



■4-2024 

arcsec 

Planetary nebula BD+30: VLA radio Image 





arcsec 



Fig. 20. Planetary nebula BD+30 imaged by the Hubble Space Telescope (top) 
and by the Very Large Array (VLA) radiotelescope in New Mexico (bottom). The 
latter indicates a complete ring of fire, dust blocks some of the visible light on the 
upper image. 

shell like structures in the time evolution of exploding fireballs comes from stellar 
astronomy. Stars with initial masses of less than about eight solar masses end their 
lives by ejecting planetary nebulae, stellar remnants turning to white dwarfs. After 
the star has completed its core hydrogen burning, it becomes a red giant. In the 
core of the star, helium burns while hydrogen continues to burn in a thin shell 
surrounding the core. This hydrogen rich shell swells to enormous size, and the 
surface temperature drops to a rather low value for stars. A solar wind develops 
that carries away most of the hydrogen envelope surrounding the star's central core. 
The envelope material ejected by the star forms an expanding shell of gas that is 
known as a planetary nebula. Planetary nebulae are illuminated by their central 
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stars and display a variety of often beautiful structures. Some are spherical or 
helical, others have bipolar shapes, and still others are rather irregularly shaped. In 
a matter of a few tens of thousands of years, they intermingle with the interstellar 
medium and disperse. 

The space-time evolution of planetary nebulae is in many aspects similar to the 
solution of non-relativistic hydrodynamics given in Ref. We argued, that this 
solution seems to describe also low and intermediate heavy ion collisions in the 30 
- 80 AMeV energy domain. A similar hydro solution may also describe the non- 
relativistic transverse dynamics at mid-rapidity in hadron -I- proton collisions in the 



CERN SPS energy domain, compare Figs |lj and the latter from Ref. [12£]. 

In all of these physical systems, expansion competes with the drop of the pres- 
sure gradients, which in turn is induced by the drop of the temperature on the 
surface. If the flow is small enough, the drop of the temperature on the surface 
results in a drop of the pressure gradients on the surface, which implies density 
pile-up. On the other hand, if the flow is strong enough, it blows away the ma- 
terial from the surface, preventing the formation of shells of fire, and an ordinary 
expanding fireball is obtained. 

Finally I note that this situation is just a special class of the more general 
solutions given in Ref. ||3^ . Arbitrary number of self-similarly expanding, simulta- 
neously existing shells of fire can be described by the general form of new class of 
exact solutions of fireball hydrodynamics [|9| . 

15. Signal of Partial Symmetry Restoration from Two- 

Pion Bose— Einstein Correlations 

In this section let me summarize Ref. |^^, where the effective intercept parameter 
of the two-pion Bose-Einstein Correlation function. A* was shown to carry a sensi- 
tive and measurable signal of partial restoration of the axial symmetry and 
the related increase of the 77' production in ultra-relativistic nuclear collisions: An 
increase in the yield of the 77' meson, proposed earlier as a signal of partial C/yi(l) 
restoration, was shown to create a "hole" in the low pt region of A* . 

In the chiral limit (m„ — ma = rris = 0), QCD possesses a U{3) chiral sym- 
metry. When broken spontaneously, U (3) implies the existence of nine massless 
Goldstone bosons. In Nature, there are only eight light pseudoscalar mesons, a 
discrepancy which is resolved by the Adler-Bell-Jackiw C/a(1) anomaly; the ninth 
would-be Goldstone boson gets a mass as a conseq uence of t he non-ze ro density of 
topological charges in the Q CD vacuum [ 141 , 14C |. In Refs | 132| , 133 1, it is argued 



that the ninth ("prodigal" |132(| ) Goldstone boson, the ij' , would be abundantly 
produced if sufficiently hot and dense hadronic matter is formed in nucleus-nucleus 



collisions. Estimates of Ref. 132 show that the corresponding production cross 
section of the rj' should be enhanced by a factor of 3 up to 50 relative to that for p 
-I- p collisions. 

If the 77' mass is decreased, a large fraction of the rj's will not be able to leave the 
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Fig. 21. Using the estimates of pion abundances given by Fritiof, the solid line 
represents A*(mt) assuming normal j]' abundances while the other lines represent 
\^{mt) with a factor of 3 (dashed), 16 (dotted) and 50 (dot-dashed) enhancement of 
rl due to partial V a(X) chiral symmetry restoration and the corresponding decrease 
of the r\' mass in the hot and dense region. All curves are calculated for Tq = 140 
MeV and (ut) = 0.5. The datapoints are from 200 AGeV central S + Pb reactions 
at CERN SPS, as measured by the NA44 Collaboration. 



0.2 
0.1 



hot and dense region through thermal fluctuation since they need to compensate for 



the missing mass by large momentum [132-134]. These 77's will thus be trapped in 



the hot and dense region until it disappears, after which their mass becomes normal 
again; as a consequence, the i]' -% will have small transverse momenta pj. Then they 
decay to pions via 

77' ?7 + 7r+ + tt" ^ (tt" + 7r+ + tt" ) + 7r+ + 7r~ . (165) 

It is important to observe that the pf of pions produced in this decay chain is small 
since many of the 77' appear at ~ and also since the rest mass of the decay 
products from the 77', decays use up most of the remaining energy. Based on the 
kinematics of the 77', ?y decay chain to pions, an enhanced production of tt mesons 
was estimated to happen dominantly in the pt ~ 150 MeV region, extending to a 
maximum pt ~ 407 MeV |6^. In the core/halo picture the vl ,t] decays contribute 
to the halo due to their large decay time (l/F,,',,, ^ 20 fm/c). Thus, we expect 



Particle Interferomctry from 40 MeV to 40 TcV 



61 



a hole in the < pt < 150 MeV region of the effective intercept parameter, A* = 

[iVeorc(p)/iVtotal(p)]'. 

To calculate the 7r+ contribution from the halo region, the bosons {uj, rj' , rj and 
Kg) are given both a rapidity (— 1.0 < ?/ < 1.0) and an nit, then are decayed using 



Jetset 7.4 [137|. The mt distribution [|l^,Q of the bosons is given by 

N{mt) = Cm^e-™*/'^'=f% (166) 

where C is a normalization constant, a = 1 — d/2 and where 

Teff = rfo + m(ui)2. (167) 

In the above expression, d = 3 is the dimension of expansion, Tfo — 140 MeV is the 
freeze-out temperature and (ut) is the average transverse flow velocity. It should be 



noted that the mt distribution for the core pions is also obtained from Eq. (166). 
The contributions from the decay products of the different regions (halo and core) 
are then added together according to their respective fractions, allowing for the 
determination of A*(m(). The respective fractions of pions are estimated separately 



by Fritiof ||l3§ and by RQMD ||l3g| as summarized in Ref. p6 |. 

Simulating the presence of the hot and dense region involves increasing the 
relative abundance of the rj' and also changing their pt spectrum. The pt spectrum 
of the rj' is obtained by assuming energy conservation and zero longitudinal motion 
at the boundary between the two phases. This conservation of transverse mass at 
the boundary implies 

m*? + Puy = + Pt^' , (168) 

where the * denotes the -q' in the hot dense region. The pt distribution then becomes 
a two-fold distribution. The first part of the distribution is from the 77' which have 
Pt ^ [^Ifi ~ rn*?]^^'^. These particles are given a pt — 0. The second part of the 
distribution comes from the rest of the ry"s which have big enough pt to leave the 
hot and dense region. These have the same, flow motivated pt distribution as the 
other produced resonances and are given a pt according to the rrit distribution 

Af^,(r7i:) = Cmr°-^e^"*/^', (169) 

where C is a normalization constant and where T' = 200 MeV and m*, is the 
effective temperature and mass, respectively, of the hot and dense region. Using 
the value given above for the effective temperature and letting m*, = 500 MeV 
implies an increase in the production cross section of the rj' in the hot and dense 
region by a factor of 10. 

Using three different effective masses for the rj' in the hot and dense region, 
calculations of A*(m() including the hot and dense regions are compared to those 
assuming the standard abundances in Fig. A similar mt dependence but with 
slightly higher values of A*(mt) is obtained when using RQMD abundances. The 
effective mass of 738 MeV corresponds to an enhancement of the production cross 
section of the rj' by a factor of 3, while the effective mass of 403 MeV and 140 MeV 
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Fig. 22. Using the estimates of pion abundances given by Fritiof, A*(rnt) is calcu- 
lated using {ut) = 0.00 (solid line), (wf) = 0.25 (dashed line), {ut) = 0.50 (dotted 
line) and {ut) = 0.75 (dashed-dotted line) 

correspond to factors of 16 and 50 respectively. The two data points shown are 
taken from NA44 data on central S + Pb reactions at the CERN SPS with incident 
beam energy of 200 AGeV ||l5| . The lowering of the -q' mass and the partial chiral 
restoration result in a hole in the effective intercept parameter at low nit- This 
happens even for a modest enhancement of a factor of 3 in the -q' production. Similar 
results are obtained when using RQMD abundances. See Ref. |^ for further details 
of the simulation. 
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In addition, A, {rrit) is calculated using Fritiof abundances with different average 
flow velocities in Fig. |2^. Here it is shown that A*(mt) can also be a measure of 
the average collective flow. In Ref. an average flow velocity of (ut) = 0.50 
resulted in an approximately flat, rrit independent shape for the effective intercept 
parameter A*(mt) distribution ||6^. Calculations using RQMD abudances result in 
a similar dependence on (ut), but with slightly higher values of A,(mf). 

This analysis of NA44 S + Pb data indicated no visible sign of ?7a(1) restoration 
at SPS energies. In addition, a mean transverse flow of (uf) w 0.50 in S + Pb 
reactions was deduced . The suggested A, -hole signal of partial restoration 
cannot be faked in a conventional thermalized hadron gas scenario, as it is not 
possible to create significant fraction of the rj and rj' mesons with ~ in such a 
case, Ref. [ p7| . 



16. Squeezed Correlations and Spectra for 
Mass-Shifted Bosons 

In this section, let me follow the lines of Refs [^[T^] to show that novel back- 
to-back correlations (BBC) arise for thermal ensembles of squeezed bosonic states 
associated with medium- modified mass shifts. It was observed in Ref. |70|, that the 
strength of the BBC could become unexpectedly large in heavy ion collisions, and 
may thus provide an experimentally observable signal of boson modification in hot 
and dense matter. 

Consider, in the rest frame of matter, the following model Hamiltonian, 

H = Ho-^J d3xd3y0(x)<5M2(x - y)(t>{y), (170) 
where Hq is the asymptotic Hamiltonian, 

The scalar field 0(x) in this Hamiltonian, H, corresponds to quasi-particles that 
propagate with a momentum-dependent medium-modified effective mass, which is 
related to the vacuum mass, mo, via 

mli\k\)=ml-SM\\k\). 

The mass shift is assumed to be limited to long wavelength collective modes: 

SM^{\k\) < ml if |k| > A,. 

The invariant single-particle and two-particle momentum distributions are given as: 

A^i(ki) = ^ki^ =c^ki (172) 
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7V2(ki,k2) = ^^kit^k2(akiaL«k2«ki>' (173) 

(«iLi "■k2"-ki) = (aici "^ki > (aL «k2 ) + (aici ^k^ > (^L «ki ) 

+ («>K«k,). (174) 

where ak is the annihilation operator for asymptotic quanta with four-momentum 
A;'' = (tJk,k), LJk = V ni^ + and the expectation value of an operator O is 
given by the density matrix p as (O) = Tr/5 0. Eq.(174) has been derived as 
a generalization of Wick's theorem for locally equilibriated (chaotic) systems in 
Rcf. p5| . 

The chaotic and squeezed amplitudes were introduced as 



Gc(l,2) = V^kiWk2(4i«k2)' 

Gs(l,2) = Vc^kit^ksiflki^ks)- 



(175) 
(176) 



In most situations, the chaotic amplitude, G'c(l,2) — 0(1,2) is dominant, and car- 
ries the Bose-Einstein correlations, while the squeezed amplitude, Gs(l, 2) vanishes. 




Fig. 23. Illustration of the new back-to-back correlations for mass-shifted 7r° pairs, 
assuming T = 140 MeV, Gc/s(pi,P2) oc exp[-Qf2-RG/2], with Rq = 2 fm. The fall 
of BBC for increasing values of |k| is controlled here by a momentum-dependent 
effective mass, TOtt,* = ^^[1 + exp(— k^/A^)], with — 325 MeV in the sudden 
approximation. 
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16.1. Mass modiGcation in a homogenous heat bath 

The terms involving Gs( 1,2) become non- negligible when mass shift becomes non- 
vanishing, i.e. dAP{]<.) 0. Given such a mass shift, the dispersion relation is 
modified to ~ "^k ~ ^M"^ (k) , where is the frequency of the in- medium mode 
with momentum k. The annihilation operator for the in-medium quasi-particle 5k, 
and that of the asymptotic field, at, are related by a Bogolyubov transformation 



1; 



where Ck = cosh[rk], Sk = sinh[rk] and rk reads as 

rk ^ log(t^k/^^k) ■ 



(177) 



(178) 



We introduce the shorthand, C-^ and Sl_i, to simplify later notation. As the Bo- 
golyubov is a squeezing transformation, let us call rk mode dependent squeeze 
parameter. While it is the a-quanta that are observed, it is the b-quanta that are 
thermalized in medium |153|. Let us consider the average for a globally thermalized 
gas of the 5-quanta, that is homogenous in volume V: 



1 



: exp 



1 V 

"t(2^ 



d^kfik blbj^ 



When this thermal average is applied, 

Gc(l,2) = ^Wki^ka 
Gs(1,2) = ^CJki^ka 



(179) 

(180) 
(181) 



If this thermal b gas freezes out suddenly at some time at temperature T, the 
observed single a-particle distribution takes the following form: 



iVi(k) 
ni(k) 
nk 



^a;kni(k). 



k -r |s_k| (n-k + 1), 
1 



exp(r!k/r) - 1 



(182) 

(183) 
(184) 



This spectrum includes a squeezed vacuum contribution in addition to the mass- 
modified thermal spectrum. 

In this homogeneous limiting case, the two particle correlation function is unity 
except for the parallel (HBT) and antiparallel (BBC) cases: 



C2(k,k) 

C2(k,-k) 



l4gk"k + C-k5-k('^-k + l)l^ 

ni(k)ni(-k) 



(185) 
(186) 
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The dynamical correlation due to the two mode squeezing associated with mass 
shifts is therefore back-to-back, as first pointed out in Ref. The strength of the 
HBT correlations remains 2 for identical momenta. 



It follows from Eq. (186) that the intercept of the BBC is unlimited from above: 
1 < C2(k, -k) < oo. As |k| oo, C2(k, -k) ~ 1 + l/|s_kP ^ 1 + l/?^l(k) ^ oo. 
Hence, at large values of |k| , the particle production is dominated by that of back-to- 
back correlated pairs for any non- vanishing value of the in-mcdium mass shifts [ [70| . 

16.2. Suppression by finite duration of emission 

To describe a more gradual freeze-out, the probability distribution F{ti) of the 
decay times ti is introduced. The sudden approximation is recovered in the F(ti) = 
6{ti — to) limiting case. The time evolution of the operators is given by ak{t) = 
o-k{ti) exp[— iLJk(^ — ti)]. This leads to a suppression of BBC as 

C2(k,-k) = l + |fK + c._0P ... 7 , T*^ — -■ (187) 

ni(k)ni(-k) 

Here F{uj) = / dtF(t) exp{—iujt), so for an exponential decay, F{t) = 8(t — to)x 
rexp[— r(i — to)] the suppression factor is 

\F{lo^ + c^_k)|' - 1/[1 + (^k + ^-k)Vr']. (188) 

In the adiabatic limit, F ^ 0, this factor suppresses completely the BBC, while 
in the sudden approximation, F — s- cx), the full strenght of the BBC is preserved. 
For a typical St — Ti/T = 2 fm/c decay time, and for BBC of (p mesons with 
TO* = 0.6 — 1.4 GeV, this suppression factor is about 0.001, which decreases the 
BBC of (j> mesons from the scale of 2000 to 2, the scale of the HBT correlations. 
This emphasizes the enormous strength of the BBC [ [70| . 

The formalism to evaluate the BBC for locally thermalized, expanding sources 
was also developed, see Ref. Q for greater details. 

As the Bogolyubov transformation always mixes particles with anti-particles, 
the above considerations hold only for particles that are their own anti-particles, 
e.g. the (j) meson and 7r°. The extension to particle-anti-particle correlations is 
straightforward. Let -I- label particles, — antiparticles if antiparticle is different 
from particle, let label both particle and antiparticle if they are identical. The 
non-trivial correlations from mass modification for pairs of (+-f), (H — ) and (00) 
type read as follows: 

r»»lk VI - I I IG.(i.2)l' , |G.(1,2)P 
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where we assume that mass modifications of particles and anti-particles are the 
same as happens at vanishing baryon density. 

This theory of particle correlations and spectra for bosons with in-medium 
mass shifts predicts huge back-to-back correlations of (j>'^,(j)'^ and K+,K^ meson 
pairs |7|]. These BBC could become observable at the STAR and PHENIX heavy 
ion experiments at RHIC [145|, and could be looked for in present CERN SPS 
experiments. Further model calculations are required to study the mass-shift effects 
on realistic source models. 



17. A Pion-Laser Model and Its Solution 

In high energy heavy ion collisions hundreds of bosons are created in the present 
CERN SPS reactions when Pb -I- Pb reactions are measured at 160 AGeV laboratory 
bombarding energy. At the RHIC accelerator, thousands of pious could be produced 
in a unit rapidity interval . If the number of pious in a unit value of phase-space 
is large enough these bosons may condense into the same quantum state and a pion 
laser could be created Q . 

In this section a consequent quantum mechanical description of multi-boson sys- 
tems is reviewed, based on properly normalized projector operators for overlapping 
multi-particle wave-packet states and a model of stimulated emission, following the 
lines of Refs [|5|,|6[|6^. One of the new analytic results is that multi-boson correla- 
tions generate momentum- dependent radius and intercept parameters even for static 
sources, as well as induce a special directional dependence of the correlation func- 
tion. This is to be contrasted to the simplistic but very frequently invoked picture 
of Eq. (H), where sources without expansion correspond to a correlation function 
that depends only on the relative momentum, but not on the mean momentum of 
the particle pairs. 

A solvable density matrix of a generic quantum mechanical system is 



/n / n \ 
A — 1 \ <„\ 1 1 / 



n=0 ^ ' " i=l \cr(") k=l 



(192) 



Here the index n characterizes sub-systems with particle number fixed to n, the mul- 
tiplicity distribution is prescribed by the set of {pnj^o' normalized as XlJ^o^*" ~ 
The density matrixes are normalized as Ti p = 1 and Trp„ = 1. The states 
\ai, an) denote properly normalized n-particle wave-packet boson states: 

I ai, - , an) = ( ^n(a.|a..) J a}, ... al\0). (193) 

Here cr^"-' denotes the set of all the permutations of the indexes {l,2,...,n} and 
the subscript sized Oi denotes the index that replaces the index z in a given per- 
mutation from cr'"'. The wave-packet creation operators, a| create the normalized 
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single-particle states \ai) = q;||0), with {ai\ai) = 1. The at = {^i,Tri,ai,ti) stands 
for a given value of the parameters of a single-particle wave-packet: the mean coor- 
dinate, the mean momentum, the width of the wave-packet in coordinate space and 
the time of the production. The distribution function pi{ai) provides the probabil- 
ity distribution for a given value of the wave-packet parameters. For simplicity, we 
assume a static source at rest, uniform wave-packet widths and simultaneous pro- 
duction, (Ti = (7 and ti = t{). A Gaussian distribution of the centers of wave-packets 
is also assumed: in the coordinate space, the distribution of is a characterized 
with a radius i?, while in the momentum-space, the centers of wave-packets tt^ are 
assumed to have a non-relativistic Boltzmann distribution corresponding to a tem- 
perature T and mass to. The coefficient of proportionality, M{n), can be determined 
from the normalization condition. 



The density matrix given in Eq. (192) describes a quantum- mechanical wave- 



packet system with induced emission, and the amount of the induced emission 
is controlled by the overlap of the n wave-packets |56| , yielding a weight in the 
range of [l,7i!]. Although it is very difficult numerically to operate with such a 
wildly fluctuating weight, the problem of overlapping multi-boson wave-packets with 
stimulated emission was reduced in Rcfs ^ to an already discovered "ring"- 
algebra of permanents for plane- wave outgoing states [^3| , with modified source 
parameters |55|,^. 

Assuming a non-relativistic, non-expanding Gaussian source at rest, and a Pois- 
son multiplicity distribution pn^ in the rare gas limiting case: 

plo) = ^exp(-no), (194) 
n! 

the ring-algebra was reduced in Ref. to a set of recurrences, which reduced the 
complexity of the problem from the numerically impossible nl to the numerically 
easy n^. These recurrences were solved analytically in Refs ||5^,Q, further reducing 
the complexity of the problem to n'', and yielding analytic insight to the behavior 
of the multi-boson symmetrization effects. 

The probability of events with fixed multiplicity n, the single-particle and the 
two-particle momentum distribution in such events are given as 



(195) 

7v(")(ki) = ^^G.(l,l), (196) 




i=l 



7v(")(ki,k2) = V V ^^[G^(l,l)Gj_„(2,2)-fG™(l,2)Gj_™(2,l)], 

(197) 



where iUn = Pn/po- Averaging over the multiplicity distribution pn yields the inclu- 
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sive spectra as 



G(l,2) - 5^G„(1,2), 

oo 

7Vi(ki) = ^p„7V(")(ki) = G(l,l), 

n=l 

Af2(ki,k2) = G(1,1)G(2,2) + G(1,2)G(2,1) 
Let us introduce the following auxiliary quantities: 
1 



7± 



[l + x±VT+2x) , 



a; = Rla'^ 



Rl = 



niT 

■1J2 ' 



(198) 

(199) 
(200) 

(201) 
(202) 



The general analytical solution of the model is given through the generating function 
of the multiplicity distribution pn 



G{z) = J2p^'" = exp £ G„(z" - 1) , 



n=0 



where G„ is introduced as 



1 



G„ = - /d^ki G„(l,l) = 



The general analytic solution for the functions G,i(l,2) is given as: 



G„(l,2) 



Jr, 



j„exp<; 



74? ki - ^1^2^ 



7+ k2 - 7- ki 



■ 7- 



(203) 

(204) 

(205) 
(206) 



The detailed proof that the analytic solution to the multi-particle wave-packet model 
is indeed given by the above equations is described in Ref. |Q . 

The representat ion of E g. (203) indicates that the quantities G„-s are the so 
called combinants |147- 14g|| of the probability distribution of p„ and in this case 
their explicit form is known for any set of model parameters. In the generator func- 
tional formalism of multi-particle production, the combinants can be introduced in 
general as the integrals of the exclusive correlation functions |15C|. The form of the 
multiphcity distribution, given by Eqs ( ^03 , ^04 ) does not correspond to the mul- 
tiplicity distributions described in standard textbooks of mathematical statistics, 
e.g. Ref. [Eol. It has the very interesting property, that the probability distribution 
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simultaneously corresponds to an infinite convolution of independently distributed 
clusters of particle singlets, pairs, triplets and higher-order n-tuples, as well as to 
an infinite convolution of strongly correlated Bose-Einstein distribution |5^, ^ of 
particle singlets, pairs, triplets etc. As far as I know, this is a new type of physically 
motivated discrete distribution in the theory of probability and statistics. 

The large n behavior of p„ depends on the ratio of no/nc, where the critical 



value of no is Uc = 7/ , ||53|-|5q]. If no < Uc, one finds (n{n — 1)) > (rt)^, a 
super-Poissonian multiplicity distribution, and a chaotic or thermal behavior of the 
inclusive correlations, C2(k, k) = 2. If no > ric, the multiplicity distribution, the 
inclusive spectra and the inclusive correlations become mathematically undefined, 
but the exlusive quantities remain finite for any fixed value of n. To calculate inclu- 
sive observables, a regularization has to be introduced similarly to the description 
of Bose-Einstein condensation of massive quanta in the limit of /i ^ m in standard 
statistical mechanics. 

Highly condensed limiting case. In Refs |^,|5^ we have related the divergence 
for no > He of the mean multiplicity (n) to the onset of a generalized type of 
Bose-Einstein condensation of the wave-packets to the wave-packet state with the 
smallest energy. Note that the onset of Bose-Einstein condensation happens in 
the limit when Pn/Pn+i ~^ 1, which happens if no Uc from below [5^ , and 
this limiting case formally corresponds to an "infinite temperature" case |151] — if 
the finite slope parameters of the 7V;["'' (k) single-particle distributions in exclusive 
events are not taken into account and the concept of the temperature is inferred 
only from the number distribution. 

In a physical situation, the total number of pions is limited: n < ue = £'tot/£'0; 
where Eq is the energy of the wave-packet with the smallest energy (including the 
mass to). Thus, energy conservation induces a cut-off in the number of pions, that 
has to be taken into account explicitly | 146| , 152 1. Such a cut in the multiplicity 
distribution can be straightforwardly implemented, as the basic building block, 
the fully symmetrized n-particle invariant momentum distribution in events with 
exactly n particles is always finite for every fixed values of n, similarly to the 
bosonic enhancement factor w„. At no — Uc, the series Sn = X]j=o changes from 
a convergent to a divergent one. After the regularization of the model, by assigning 
a zero probability to multiplicities greater that Ue, one can show that for no > nc a 
Bose-Einstein condensation develops more and more with increasing values of no. 

Utmost care is required when evaluating the results published in the literature 
regarding the nature of coherence and Bose-Einstein condensation in the pion-laser 
model: some papers identify the "Bose-Einstein condensation" with the "infinitely 
hot" no — He limiting case. At this point, however, the condensate just appears 
with non-zero probability (and one has to introduce the cut multiplicity distribution 
to describe it with a > 0), but the number of quanta in the condensate is rather 
small, PriE c< 1/n^; at the n = ric critical point. 

The nature of the Bose-Einstein condensation was discussed and clarified in 
Ref . |152| , where it was shown that the condensate will fully develop and dominate 
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the density matrix in the R ^ and T — > simuhaneous hmiting cases, confirming 
the intuitive picture that Bose-Einstein condensation happens in very cold and very 
small systems. 

In the highly condensed limiting case, the multiplicity distribution of the pro- 
duced particles will be sub-Poissonian, a very narrow, cut power-law distribution 
that increases with n as 

Pn oc 1 — ) 0(n_E - n) for Uq > n^, (207) 

and vanishes after n > ue- In the limit when the number of particles in the 
condensate is very large, the exclusive and the inclusive correlation functions become 
unity |§, 

C(ki,k2) = C("nki,k2) = l (208) 
(the highly condensed limiting case, hq 3> ric). 

By definition, the above equalities imply optical coherence in the highly condensed 
limiting case ^6, 152|. It is worthwhile to emphasize, that optical coherence is 



not to be confused with the appearance of the coherent states of the annihilation 



operator |152|. Instead of being an eigenstate of the annihilation operator, the fully 
developed Bose-Einstein condensate is an eigenstate of the creation operator, with 
zero eigenvalue. This is due the cutoff induced by the conservation of energy: it is 
not possible to add one more pion to the condensate if already all the pions allowed 
by the constraints are in the condensate. 

Rare gas limiting case. In contrast, the large source sizes or large effective 
temperatures correspond to a rare Boltzmann gas, the a; 3> 1 limiting case. The 
general analytical solution of the model becomes particularly simple in this limiting 
case. The leading order multiplicity distribution can be found from Eqs ( p03| , p04| ) 



corresponding to independently distributed particles with a small admixture of in- 
dependently distributed particle pairs |5q]: 



Pn = — |-exp(-no) 
n! 



n{n — 1) - 
2(2a;)^ 



(209) 



The mean multiplicity, the factorial cumulant moments of the multiplicity distribu- 
tion, the inclusive and exclusive momentum distributions were obtained to leading 
order terms in 1/x in Rcf. |^5|] . Figure ^ indicates that the radius parameter of 
the exclusive correlation function becomes mean momentum momentum- dependent, 
even for static sourcesl This genuine multi-particle symmetrization effect is more 
pronounced for higher values of the fixed multiplicity n, in contrast to the momen- 
tum dependence of Ak that is independent of n [^5|| . 

One finds that multi-boson symmetrization effects lead to the development of 
a Bose-Einstein condensate. Before the onset of the Bose-Einstein condensation, 
the stimulated emission becomes significant in the low momentum modes earlier 
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Multi-Particle Symmetrization Effects 
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Fig. 24. Multi-particle symmetrization results at low K in a momentum-dependent 
reduction of the intercept parameter Ak, the side- wards and the outwards radius 

parameters, Rk.,s and Rk.o from their static values of 1 and Re, respectively. The 
enhancement of these parameters at high momentum is hardly noticeable for large 
and hot systems. 

than in the high momentum modes. This is the reason, why even the exclusive 
correlation functions develop a mean momentum dependent radius parameter, as 
well as a direction-dependent radius component and a mean momentum dependent 
intercept parameter. 

18. Summciry and Outlook 

In this review, new kind of similarities were highlighted between stellar astronomy 
and intensity interferometry in high energy physics. The model independent charac- 
terization of short-range correlation was given in terms of expansions in complete or- 
thonormal sets of polinomials, the core/halo model and the recently found Coulomb 
wave-function correction method was reviewed, for Bose-Einstein n-particle corre- 
lations. 

The invariant Buda-Lund (BL) parameterization of Bose— Einstein correlation 
functions was derived in a general form, and compared to the Bertsch-Pratt and 
the Yano-Koonin-Podgoretskii parameterization in the particular Gaussian limit- 
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ing case. The Buda-Lund hydrodynamical parameterization, BL-H was fitted to 
hadron-proton and Pb + Pb collisions at CERN SPS energies. Larger mean freeze- 
out proper-times and larger transverse radii were found in the Pb -I- Pb reactions. 
Although the central values of freeze-out temperatures were rather similar in both 
reactions, the transverse temperature gradient is larger while the transversal flow is 
smaller in h -I- p reactions, than in the Pb + Pb system. This resulted in different 
shapes for the transverse density profiles, that were approximately reconstructed as- 
suming the applicability of a new family of solutions to fireball hydrodynamics [ p9[ . 
Although Pb -|- Pb reactions were found to be rather homogenous expanding fire- 
balls, the h -|- p reactions were found to be similar to a cold and expanding ring 
of fire when viewed in the transverse plane. The central freeze-out temperature is 
about To = 140 MeV in both reactions, the surface temperature after the emission 
of particles is over seems to be also similar, about Tg = 82 MeV, the duration of 
the particle emission is also about Ar « 1.5 fm in both cases. 

Inspecting the results of a non-relativistic version of BL-H to '"'Ar 4- ^^^Au 
proton and neutron correlations and spectra, an indirect signal was observed for 
the formation of a shell of fire, made of protons, while the neutrons seems to come 
from an ordinary fireball. The hydrodynamics of cooling and expanding shells of 
low energy heavy ion reactions was shown to be similar to that of spherical plan- 
etary nebulae, indicating a new connection between stellar astronomy and particle 
interferometry in heavy ion physics. 

Another similarity between stellar astronomy and high energy physics was dis- 
cussed in terms of the interferometry of binary sources: the binary stars in stellar 
astronomy create oscillations in the HBT effect [113| similarly to the oscillations 
that were shown to exist in the Buda-Lund type of hydrodynamical parameteri- 
zation in heavy ion physics and to the expected oscillations of pion correlations in 
particle interferometry in W+W^ decays at LEP2. The first positive evidence for 
the existence of such binary sources in heavy ion physics seems to be the recent mea- 
surement of oscillating proton-proton correlations by the NA49 Collaboration |114[ , 
which may be a consequence of the existence of the two maxima in the proton rapid- 
ity distribution and the attractive final-state interactions of protons, that enhance 
the large Q part of the pp intensity correlation function and make these oscillations 
clearly visible. 

The question of non-Gaussian oscillations of three-dimensional Bose-Einstein 
correlation functions in heavy ion physics has not yet been experimentally investi- 
gated. I think it is time to start the experimental search for non-Gaussian structures 
in multi-dimensional Bose-Einstein correlation functions in high energy heavy ion 
and particle physics. I hope that experiments will decide to publish in the fu- 
ture not only the (Gaussian) fit parameters of (multi-dimensional) Bose-Einstein 
correlation functions, but, most importantly, the measured data points and the 
corresponding the error bars. It was shown already in Refs |l^-[20|] , that the recon- 
struction of the space-time picture of the particle emission: the extraction of den- 
sity, fiow and temperature profiles requires the simultaneous analysis of the double- 
differential single-particle spectra and the momentum-dependent multi-dimensional 
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Bose-Einstein correlation functions. These data sets should be made pubUc in as 
much detail as possible, including multi-dimensional data and error-bar tables. 

In the chiral limit, when the up, down and strange quarks become massless 
and the UaI^) symmetry is fully restored, the mass of the rj' meson vanishes in 
the ^/^(l) symmetric, new phase. The appearance of such a phase implies that 
the intercept parameter of the two-pion correlation function vanishes in the pt < 
150 MeV region. In this sense, the transverse mass-dependent intercept parameter 
A*(mf), was interpreted as an effective order parameter of partial [/^(l) symmetry 
restoration p7 , 



Bosonic mass shifts in medium were shown to result in unlimitedly large back- 
to-back correlations of the observable boson-anti-boson pairs. Although a finite 
time suppression factor may reduce the strength of these correlations substantially, 
the magnitude of the back-to-back correlations is estimated to be observably strong 
for typical mass shifts and freeze-out time distributions in ultra-relativistic heavy 
ion collisions. 

Multi-boson symmetrization effects were shown to generate momentum- depend- 
ent radius and intercept parameters even for static sources. 

The proposed A*-hole signal of the Ua{^) symmetry restoration, the new kind 
of back-to-back correlations and optically coherent, effectively lasing pion sources 
could be searched for in future in heavy ion experiments at CERN SPS and at RHIC. 
The oscillations in multi-dimensional Bose-Einstein and Fermi-Dirac correlations 
could be searched for in e+e^ annihilation experiments at LEP2, as well as in heavy 
ion collisions at CERN SPS and at RHIC. Note that this paper is a review of particle 
interferometry before 2000; a substantially shortened version of the present material 
has been published in Ref. |154|. 
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